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ABSTRACT 

We present echelle spectrophotometry of the blue compact dwarf galaxy (BCDG) NGC 5253. 
The data have been taken with the Very Large Telescope UVES echelle spectrograph in the 3100 
to 10400 A range. We have measured the intensities of a large number of permitted and forbidden 
emission lines in four zones of the central part of the galaxy. In particular, we detect faint C II 
and O II recombination lines. This is the first time that these lines are unambiguously detected 
in a dwarf starburst galaxy. The physical conditions of the ionized gas have been derived using 
a large number of different line intensity ratios. Chemical abundances of He, N, O, Ne, S, CI, 
Ar, and Fe have been determined following the standard methods. In addition, C"^"*" and 0+"'" 
abundances have been derived from pure recombination lines. These abundances are larger than 
those obtained from collisionally excited lines (from 0.30 to 0.40 dex for C''"+ and from 0.19 
to 0.28 dex for 0+"'"). This result is consistent with a temperature fluctuations parameter (t^) 
between 0.050 and 0.072. We confirm previous results that indicate the presence of a localized 
N enrichment in certain zones of the center of the galaxy. Moreover, our results also indicate a 
possible slight He overabundance in the same zones. The enrichment pattern agrees with that 
expected for the pollution by the ejecta of massive stars in the Wolf-Rayet (WR) phase. The 
amount of enriched material needed to produce the observed overabundance is consistent with 
the mass lost by the number of WR stars estimated in the starbursts. Finally, we discuss the 
possible origin of the difference between abundances derived from recombination and collisionally 
excited lines (the so-called abundance discrepancy problem) in H II regions, finding that a recent 
hypothesis based on the delayed enrichment by SNe ejecta inclusions seems not to explain the 
observed features. 
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1. Introduction 

Deep cross-dispersed echelle spectrophotometry 
with large aperture telescopes is a powerful tech- 
nique for refining our knowledge about the che- 
mical composition of H II regions. It permits to 
observe the whole optical spectral range, to mea- 
sure crucial faint auroral and recombination lines, 
to deblend nebular lines, and to decontaminate 
them for nearby sky spectral features. 

Our group has obtained deep high resolution 
spectra of most of the brightest Galactic H II re- 
gions (Esteban et al. 1998, 1999a,b, 2004; Garcfa- 
Rojas et al. 2004, 2005, 2006) that have led to the 
precise determination of the abundance of heavy 
element ions from recombination lines and even 
to the determination of the carbon radial abun- 
dance gradient in the Galactic disk (Esteban et 
al. 2005). There is still a limited number of simi- 
lar studies devoted to giant extragalactic H II re- 
gions (hereafter GEHRs), focusing the efforts on 
the Magellanic Clouds (Peimbert 2003; Tsamis et 
al. 2003) and in a few other nearby irregular and 
spiral galaxies (Esteban et al. 2002; Peimbert et 
al. 2005). A common result of both the Galactic 
and extragalactic works, is the finding of a sys- 
tematic difference between the abundances deter- 
mined from recombination lines (hereafter RLs) 
and coUisionally excited lines (hereafter CELs) of 
the same ions, in the sense that abundances de- 
termined from RLs are always higher than those 
determined from CELs. This problem, that has 
been also reported to occur in some planetary ne- 
bulae (PNe), is known as the abundance discre- 
pancy. The origin of this problem is still un- 
known. Pequignot et al. (2002) and Tsamis et 
al. (2004) have constructed models with inhomo- 
geneous chemical composition and physical condi- 
tions to explain that discrepancy in PNe. A simi- 
lar idea has been very recently proposed by Tsamis 
& Pequignot (2005) for explaining the abundance 
discrepancy in the case of 30 Dor. These authors 
propose the presence of metal enhanced inclusions 
that would be responsible for most of the emis- 
sion in RLs of heavy element ions. However, the 
abundance discrepancy could be related to an- 
other secular problem in nebular astrophysics: the 
existence of temperature fluctuations in the ion- 
Based on observations collected at the European Southern 
Observatory, Chile, proposal number ESO 70.C-0008(A) 



ized gas volume (Peimbert 1967; Torres-Peimbert, 
Peimber & Daltabuit 1980), whose existence is 
still controversial although some mechanisms have 
been proposed (see Esteban 2002 for a review). 

One step ahead in our investigation is to obtain 
abundance determinations from RLs in a sample 
of bright GEHRs in starburst galaxies and even 
H II galaxies of chemical compositions other than 
solar. The dwarf irregular galaxy NGC 5253 is 
one of the closest starbursts and a very appro- 
priate target for our purposes (Figure 1). It lies 
at a heliocentric distance of 3.3 Mpc (Gibson et al. 
2000) and belongs to the Centaurus Group. Con- 
sequently, NGC 5253 has been one of the most 
studied starbursts, observed at basically all wave- 
lengths from radio to X-rays. 

NGC 5253 is also especially interesting because 
it is the best candidate for localized chemical en- 
richment in a GEHR. Welch (1970), Walsh & Roy 
(1987, 1989), and Kobulnicky et al. (1997) re- 
ported the presence of a strong nitrogen overabun- 
dance in a particular zone of the starbursting nu- 
cleus of the galaxy. Campbell et al. (1986) also 
found a helium enhancement in that zone, but this 
has not been confirmed in later works. Another 
peculiarity of the galaxy is that the majority of 
the gas seems to rotate about the optical major 
axis of the galaxy (Kobulnicky & Skillman 1995) 
but this behaviour is not clear and could be com- 
bined with some kind of outflow (Koribalski 2006, 
priv. communication). However, CO observations 
by Turner, Beck, & Hurt (1997) and Meier, Turner 
& Beck (2002) suggest that molecular clouds are 
infalling into NGC 5253. These two observational 
facts indicate that the dynamical situation of this 
galaxy is far from clear. Other authors as van den 
Bergh (1980) and CardweU & Phillips (1989) have 
suggested that the nuclear starburst was triggered 
as a result of a past interaction (around 1 Gyr ago) 
with the neighboring galaxy M83. Both galaxies 
lie at a radial distance of only 500 kpc and the pro- 
jected separation is 130 kpc (Thim et al. 2003). 

Concerning the kinematics of the ionized gas of 
NGC 5253, it has to be highlighted the presence of 
two large superbubbles related to the central re- 
gion of the galaxy, with diameters of the order of 
1 arcmin (~1 kpc), and expansion velocities of 35 
kms~^ (Marlowe et al. 1995). These superbubbles 
coincide with extended X-ray emission as demon- 
strated by Chandra and XMM-Newton data ana- 



2 



NGC 5253 - Starburst Galaxy 




Fig. 1. — (Left) False color image in filters B (blue), V (green) and / (red) of the starburst galaxy NGC 5253 observed 
by the 2.5m du Pont telescope, Las Campanas Observatory. (Right) False color image obtained by combining HST 
Ha (red), [O III] (green) and UV (blue) images. The position of our VLT slit is indicated in it. North is up and east 
at left. The HST image has been included at the center of the 2.5m du Pont telescope image for comparison. 



lyzed by Summers et al. (2004). The structure and 
morphology of the ionized gas in NGC 5253 were 
deeply analyzed by Calzetti et al. (2004) using im- 
agery from Hubble Space Telescope (HST), detect- 
ing faint arches and filaments in both Ha and [S II] 
at >1 kpc from the main ionizing cluster (see Fig- 
ure 1, right panel) that are partially excited by 
shocks. 

NGC 5253 has been considered as one of the 
youngest starbursts in the local universe (van 
den Bergh 1980; Moorwood & Glass 1982; Rieke 
et al. 1988). CampbeU, Terlevich, & Melnick 
(1986) and Walsh & Roy (1987) reported a broad 
emission feature, indicating the presence of Wolf- 
Rayet (WR) stars in this galaxy. Several authors 
(Schaerer et al. 1997; Kobulnicky et al. 1997) have 
confirmed the presence of late-type WN and early- 
type WC stars. This fact, and the almost en- 
tirely thermal radio spectrum with very little syn- 
chrotron emission from supernova remants (Beck 
et al. 1996), imply the extreme youth of the star- 
burst. Indeed, there is evidence in the literature 
that multiple starbursts in succession have oc- 
curred in the center of this galaxy, leaving behind 
an evolved, but still young-ish stellar population 
(Strickland & Stevens 2000; Tremonti et al. 2001). 



2. Observations and data reduction 

The observations were made on 2003 March 30 
with the Ultraviolet Visual Echelle Spectrograph, 
UVES (D'Odorico et al. 2000), at the Very Large 
Telescope, VLT, Kueyen unit in Cerro Paranal Ob- 
servatory (Chile). The standard settings in both 
the red and blue arms of the spectrograph, cover- 
ing the region from 3100 to 10400 A, were used. 
The wavelength intervals 5783-5830 A and 8540- 
8650 A were not observed due to a gap between 
the two CCDs used in the red arm. There are also 
five small intervals that were not observed, 9608- 
9612 A, 9761-9767 A, 9918-9927 A, 10080-10093 
A and 10249-10264 A, because the five redmost 
orders did not fit completely within the CCD. The 
whole spectrum was taken in two different blocks 
of observations. The first block covers 3800-5000 
A with the blue arm and 6700-10400 A with the 
red arm and was taken in three consecutive 1000s 
exposures. The second block, covering 3100-3900 
A and 4750-6800 A, was taken in two consecutive 
360s exposures. None of the emission lines was 
saturated in the exposures. The journal of the 
observations can be found in Table 1. 

The atmospheric dispersion corrector (ADC) 
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Table 1: Journal of observations. 



AA 


Exp. time 


Spectral res.* 


Spatial res 


(A) 


(s) 


(A pix-i) 


(" Pix-^) 


3100-3900 


2 X 360 


0.019 


0.25 


3800-5000 


3 X 1000 


0.022 


0.25 


4750-6800 


2 X 360 


0.024 


0.18 


6700-10400 


3 X 1000 


0.033 


0.17 



^ At the center of each wavelength range. 



was \iscd to keep the same observed region within 
the slit regardless of the air mass value. This is 
especially important for this work because we ex- 
tracted, analyzed and compared different small ar- 
eas along the slit. The slit width was set to 1.5" 
and the slit length was set to 10" in the blue arm 
and to 12" in the red arm. The slit width was 
chosen to maximize the S/N ratio of the emission 
lines, to separate most of the relevant faint lines, 
and to obtain a good spectral resolution for ana- 
lyzing the velocity structure of the ionized gas. 
The effective resolution at a given wavelength is 
approximately A A ^ A/ 17600. The seeing was ex- 
cellent during the observations, ~0.5 ". 

In Figure 2 we show our slit position over the 
Figure 1 of Kobulnicky et al. (1997). The slit po- 
sition was located along the north-south direction 
(P. A. = 0°), which was chosen in order to observe 
the most interesting regions inside the starburst. 
These zones were previously analyzed by Walsh & 
Roy (1989) and Kobulnicky et al. (1997). From 
north to south, these regions were designated HII- 
2, HII-1, UVl, and UV2 (A, B, C, and D, respec- 
tively), as it is shown in Figure 2. We have ex- 
tracted 1-D spectra of the regions, with a size of 
1.5" X 1.5". 

The spectra were reduced using the IRAF^ 
echelle reduction package, following the standard 
procedure of bias subtraction, flatfielding, wave- 
length calibration, aperture extraction and flux 
calibration. The standard stars EG 247, HD 
49798, and C-32d9927 were observed for flux cali- 
bration. As an example, we show the wavelength- 
and flux-calibrated spectra of region B (HII-1) in 
Figure 3. 

^IRAF is distributed by NOAO which is operated by AURA 
Inc., under cooperative agreement with NSF 
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Fig. 2. — Our VLT slit position over Figure 1 of Ko- 
bulnicky ct al. (1997), that shows the Ha contours 
over the 2200 A continuum image (gray scale). The 
four analyzed knots inside NGC 5253 are indicated 
and labeled. 

3. Line intensities 

and reddening correction 

Line intensities and equivalent widths were 

measured integrating all the flux in the line be- 
tween two given limits and over a local continuum 
estimated by eye. In the cases of line blending, 
a mtiltiple Gaussian profile fit procedure was ap- 
plied to obtain the line flux of each individual 
line. The measurements were performed with the 
SPLOT routine of the IRAF package. 

In Table 2 we have compiled all the emission 
lines detected in the four analyzed regions. Knot 
B (HII-1) is the zone where more lines have been 
identified: 169. Therefore, NGC 5253 is nowa- 
days the dwarf starburst galaxy where more opti- 
cal emission lines have been reported. We identify 
156 fines in A (HII-2), 161 in C (UVl), and 86 in D 
(UV2). Besides the H 1 Paschen lines included in 
Table 2, some other lines were also detected, but 
they are blended with telluric spectral features, 
making it impossible to obtain a good determi- 
nation of their intensity. The identification and 
adopted laboratory wavelength of the lines, as well 
as their errors, were obtained following Garcia- 
Rojas et al. (2004); Esteban et al. (2004). Colons 
indicate errors of the order or greater than 40%. 
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Table 2 

Dereddened line intensity ratios with respect to I(H,'?) = 100. 





on 


Mult 


/(A) 


(A) HII-2 


(B) HII-1 


) U V -1 




3187.84 


Ho I 


3F 





308 


4.01 ± 0.33 


2.89 ± 0.28 


3.82 ± 0.26 




3447.49 


Ho I 


7 





295 






0.309 








3450.39 


[Fo II] 


27F 





295 






0.063 








3530.50 


Ho I 


36 





286 






0.152 








3554.42 


Ho I 


34 





283 


0.249 




0.158 




0.44 ± 0.11 




3587.28 


He I 


32 





278 


0.322 




0.191 








3613.64 


He I 


6 





275 


0.335 




0.356 




0.34 ± 0.10 




3634.25 


He I 


28 





272 


0.313 




0.411 




0.28 ± 0.10 




3651.97 


Ho I 


27 





269 






0.093 








3664.68 


H I 


H28 





267 










0.24 ± 0.09 




3666.10 


H I 


H27 





267 


0.159 




0.124 




0.35 ± 0.10 




3667.28 


H I 


H26 





266 


0.130 




0.144 




0.276 ± 0.097 




3669.47 


H I 


H25 





266 


0.172 




0.183 




0.271 ± 0.097 




3671.48 


H I 


H24 





266 


0.166 




0.205 




0.256 ± 0.095 




3673.76 


H I 


H23 





265 


0.187 




0.255 




0.51 ± 0.12 




3676.37 


H I 


H22 





265 


0.284 




0.316 




0.51 ± 0.12 




3679.36 


H I 


H21 





265 


0.374 




0.229 




0.57 ± 0.12 




3682.81 


H I 


H20 





264 


0.365 




0.389 




0.68 ± 0.13 




3686.83 


H I 


H19 





263 


0.68 ± 0.20 


0.56 ± 0.18 


0.71 ± 0.13 




3691.56 


H I 


H18 





263 


0.91 ± 0.22 


0.76 ± 0.19 


0.97 ± 0.15 


0.570: 


3697.15 


H I 


H17 





262 


1.18 ± 0.23 


1.01 ± 0.21 


1.21 ± 0.16 


1.021: 


3703.86 


H I 


H16 





260 


1.30 ± 0.23 


1.26 ± 0.22 


1.24 ± 0.16 


0.736: 


3705.04 


He I 


25 





260 


0.58 ± 0.19 


0.414: 


0.49 ± 0.12 




3711.97 


H I 


H15 





259 


1.40 ± 0.24 


1.27 ±0.22 


1.41 ± 0.17 


1.324: 


3721.83 


[s ni] 


2F 





257 


2.95 ± 0.29 


3.12 ± 0.28 


3.24 ± 0.24 


2.96 ± 0.69 


3726.03 


[O II] 


IF 





257 


65.3 ± 2.3 


62.9 ± 2.2 


80.3 ± 2.9 


145.4 ± 5.3 


3728.82 


[O II] 


IF 





256 


71.9 ± 2.6 


71.9 ± 2.6 


102.2 ± 3.6 


208.0 ± 7.2 


3734.17 


H I 


H13 





255 


2.62 ± 0.28 


2.16 ± 0.25 


2.31 ± 0.20 


2.60 ± 0.67 


3750.15 


H I 


H12 





253 


3.14 ± 0.30 


2.89 ± 0.27 


2.80 ± 0.22 


2.78 ± 0.68 


3770.63 


H I 


Hll 





249 


3.96 ± 0.32 


3.62 ± 0.29 


3.90 ± 0.26 


4.59 ± 0.80 


3797.90 


H I 


HIO 





244 


5.32 ± 0.36 


5.08 ± 0.33 


5.13 ± 0.30 


5.33 ± 0.84 


3819.61 


He I 


22 





240 


1.10 ± 0.22 


0.93 ± 0.20 


0.94 ± 0.14 




3835.39 


H I 


H9 





237 


7.39 ± 0.28 


7.37 ± 0.27 


7.44 ± 0.28 


7.27 ± 0.48 


3856.02 


Si II 


1 





233 


0.079 ± 0.029 


0.204 ± 0.031 


0.086 ± 0.031 




3862.59 


Si II 


1 





232 


0.073 ± 0.029 


0.084 ± 0.022 


0.109 ± 0.034 


0.52 ± 0.17 


3868.75 


[No HI] 


IF 





230 


48.8 ± 1.7 


47.9 ± 


1.6 


34.1 ± 1.2 


24.5 ± 1.0 


3871.82 


Ho I 


60 





230 






0.092 ± 0.023 


0.038: 




3889.05 


H I 


H8 





226 


19.19 ± 0.67 


17.91 ± 0.63 


18.15 ± 0.64 


18.67 ± 0.86 


3926.53 


He I 


58 





219 






0.070 ± 0.021 


0.135 ± 0.036 




3964.73 


He I 


5 





211 


0.714 ± 0.063 


0.680 ± 0.052 


0.637 ±0.063 


0.336: 


3967.46 


[Ne HI] 


IF 





210 


13.81 ± 0.48 


14.40 ± 0.50 


9.51 ± 0.35 


6.24 ± 0.46 


3970.07 


H I 


H7 





210 


15.98 ± 0.55 


16.03 ± 0.55 


16.10 ± 0.56 


15.83 ± 0.75 


4009.22 


He I 


55 





202 


0.134 ± 0.034 


0.232 ± 0.032 


0.130 ± 0.035 




4026.21 


He I 


18 





198 


1.93 ± 0.10 


1.815 ± 0.090 


1.84 ± 0.10 


1.51 ± 0.25 


4068.60 


[S II] 


IF 





189 


2.37 ± 0.12 


1.602 ± 0.082 


1.484 ± 0.091 


3.46 ± 0.35 


4069.62 


O II 


10 





189 










0.054: 




4072.15 


O II 


10 





188 










0.020: 




4076.35 


[S II] 


IF 





187 


0.757 ± 0.064 


0.490 ± 0.044 


0.471 ± 0.055 


1.17 ± 0.23 


4097.26 


O II 


20-48 





183 


0.031 












4101.74 


H I 


H6 





182 


26.00 ± 0.86 


26.20 ± 0.86 


26.05 ± 0.86 


25.9 ± 1.0 


4120.82 


He I 


16 





177 


0.152 ± 0.036 


0.188 ± 0.030 


0.094 ± 0.031 




4143.76 


He I 


53 





172 


0.212 ± 0.040 


0.209 ± 0.031 


0.199 ± 0.040 




4168.97 


He I 


52 





167 






0.034 




0.039: 




4243.97 


[Fo II] 


21F 





149 


0.068 ± 0.027 


0.069 ± 0.020 


0.037: 


0.124: 


4267.15 


C II 


6 





144 


0.063 ± 0.026 


0.074 




0.062 ± 0.026 




4276.83 


[Fo II] 


21F 





142 


0.023 




0.029 




0.019: 




4287.40 


[Fe II] 


7F 





139 


0.143 ± 0.034 


0.108 ± 0.023 


0.184 ± 0.039 


0.40 ±0.15 


4303.61 


O II 


66 





135 










0.025: 




4340.47 


H I 


H7 





127 


47.0 ± 1.5 


46.8 ± 1.5 


46.9 ± 1.5 


46.9 ± 1.7 
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Table 2 — Continued 



Ao (A) 


Ion 


Mult. 


/(A) 


(A) HII-2 


(B) HII-1 


(C) UV-1 


(D) UV-2 


4359.34 


[Fo II] 


7F 


0.122 


0.067: 


0.103 ± 0.023 


0.099 ± 0.031 


0.279: 


4363.21 


[O III] 


2F 


0.121 


6.46 ± 0.23 


6.70 ± 0.23 


3.95 ± 0.16 


2.61 ± 0.31 


4368.25 


O I 


5 


0.120 


0.019: 


0.068 ± 0.019 


0.037: 




4387.93 


He I 


51 


0.115 


0.419 ± 0.049 


0.474 ± 0.042 


0.403 ± 0.050 


0.322: 


4413.78 


[Fe II] 


7F 


0.109 




0.095 ± 0.022 


0.073 ± 0.028 




4416.27 


[Fe II] 


6F 


0.109 


0.052: 


0.045 ± 0.017 


0.043: 




4437.55 


He I 


50 


0.104 


0.052: 


0.044 ± 0.016 


0.046: 




4452.11 


[Fe II] 


7F 


0.100 


0.049: 


4.20 ± 0.16 


0.032: 




4471.48 


He I 


14 


0.095 


4.11 ± 0.16 


4.09 ± 0.15 


3.79 ± 0.15 


3.36 ± 0.34 


4562.60 


Mg I] 


1 


0.073 


0.152 ± 0.034 


0.108± 0.022 


0.129 ± 0.033 


0.49 ± 0.16 


4571.20 


Mg I] 


1 


0.071 


0.139 ± 0.033 


0.102± 0.022 


0.115 ± 0.032 


0.290: 


4638.86 


O II 


1 


0.055 


0.040: 




0.049: 




4641.81 


O II 


1 


0.054 


0.066: 




0.060 ± 0.019 




4649.13 


O II 


1 


0.052 


0.063 ± 0.020 


0.063 ± 0.021 


0.062 ± 0.019 




4650.84 


O II 


1 


0.052 


0.034: 


0.036: 


0.042: 




4658.10 


[Fe III] 


3F 


0.050 


1.185 ± 0.074 


1.134 ± 0.062 


0.815± 0.065 


2.11 ± 0.28 


4661.63 


O II 


1 


0.049 


0.027: 


0.026: 


0.061: 




4701.53 


[Fo III] 


3F 


0.039 


0.291 ± 0.042 


0.339 ± 0.035 


0.195 ± 0.038 


0.340: 


4711.37 


[Ar IV] 


IF 


0.037 


0.681 ± 0.058 


0.890 ± 0.054 


0.170 ± 0.036 




4713.14 


He I 


12 


0.037 


0.578 ± 0.054 


0.587 ± 0.044 


0.396 ± 0.049 




4733.93 


[Fe III] 


3F 


0.031 


0.090 ± 0.029 


0.093 ± 0.021 


0.048: 




4740.16 


[Ar IV] 


IF 


0.030 


0.723 ± 0.059 


0.918 ± 0.055 


0.135 ± 0.033 




4754.83 


[Fe III] 


3F 


0.026 


0.176 ± 0.035 


0.157 ± 0.025 


0.102 ± 0.030 


0.240: 


4769.60 


[Fe III] 


3F 


0.023 


0.131 ± 0.032 


0.072 ± 0.019 


0.058: 




4814.55 


[Fe II] 


20F 


0.012 


0.055: 


0.044 ± 0.016 


0.039: 




4861.33 


H I 


H/3 


0.000 


100.0 ± 3.0 


100.0 ± 3.0 


100.0 ± 3.0 


100.0 ± 3.2 


4881.00 


[Fo III] 


2F 


-0.005 


0.345 ± 0.044 


0.291 ± 0.032 


0.170 ± 0.036 


0.340: 


4889.70 


[Fo II] 


3F 


-0.007 


0.028: 








4921.93 


Ho I 


48 


-0.015 


1.047 ± 0.068 


1.016 ± 0.057 


0.980 ± 0.069 


0.82 ± 0.19 


4931.32 


[O III] 


IF 


-0.017 


0.083 ± 0.027 


0.056 ± 0.017 


0.030: 




4958.91 


[O III] 


IF 


-0.024 


204.0 ± 6.1 


206.2 ± 6.2 


160.9 ± 4.8 


104.3 ± 3.3 


4985.90 


[Fe III] 


2F 


-0.031 


0.436 it 0.061 


0.459 ± 0.055 


0.503 ± 0.051 


1.85 zb 0.29 


5006.84 


[O III] 


IF 


-0.036 


579 ± 17 


597 ± 18 


460 ± 13 


300.2 ± 9.3 


5015.68 


Ho I 


4F 


-0.038 


1.95 ± 0.12 


2.18 ± 0.12 


1.98 ± 0.11 


1.76 ± 0.29 


5041.03 


Si II 


5 


-0.044 


0.144 ± 0.039 


0.205 ± 0.039 


0.124 ± 0.026 


0.53 ± 0.18 


5047.74 


Ho I 


47 


-0.046 


0.277 ± 0.051 


0.327 ± 0.048 


0.311 ± 0.040 


2.24 ± 0.32 


5055.98 


Si II 


5 


-0.048 


0.346 ± 0.055 


0.223 ± 0.040 


0.300 ± 0.039 




5158.81 


[Fe II] 


19F 


-0.073 




0.118 ± 0.031 


0.181 ± 0.031 


0.382: 


5191.82 


[Ar III] 


3F 


-0.081 


0.158 ± 0.040 


0.153 ± 0.034 


0.089 ± 0.022 


0.297: 


5197.90 


[N I] 


IF 


-0.082 


0.318 ± 0.053 


0.301 ± 0.045 


0.220 ± 0.034 


0.76 ± 0.20 


5200.26 


[N I] 


IF 


-0.083 


0.251 ± 0.048 


0.237 ± 0.041 


0.147 ± 0.028 




5261.61 


[Fe II] 


19F 


-0.098 


0.074: 


0.042: 


0.061 ± 0.019 




5270.40 


[Fe III] 


IF 


-0.100 


0.466 ± 0.061 


0.445 ± 0.053 


0.324 ± 0.040 


0.77 ± 0.21 


5517.71 


[CI III] 


IF 


-0.154 


0.375 ± 0.055 


0.320 ± 0.045 


0.362 ± 0.042 


0.254: 


5537.88 


[CI III] 


IF 


-0.158 


0.289 ± 0.050 


0.251 ± 0.040 


0.271 ± 0.037 


0.105: 


5754.64 


[N II] 


3F 


-0.194 


0.500 ± 0.062 


0.439 ± 0.050 


0.177 ± 0.030 


0.316: 


5875.64 


He I 


11 


-0.215 


11.50 ± 0.43 


12.34 ±0.45 


11.14 ± 0.42 


10.74 ± 0.65 


5957.56 


Si II 


4 


-0.228 


0.160 ± 0.039 








5978.93 


Si II 


4 


-0.231 




0.035: 


0.066 ± 0.019 




6300.30 


[OI] 


IF 


-0.282 


2.32 ± 0.13 


2.15 ± 0.12 


2.35 ± 0.12 


6.05 ± 0.48 


6312.10 


[S III] 


3F 


-0.283 


2.51 ± 0.14 


2.43 ± 0.13 


1.705 ± 0.099 


1.53 ± 0.26 


6347.11 


Si II 


2 


-0.289 


0.119 ± 0.034 


0.076 ± 0.023 


0.108 ± 0.023 




6363.78 


[O I] 


IF 


-0.291 


0.730 ± 0.072 


0.677 ± 0.060 


0.739 ± 0.060 


1.90 ± 0.29 


6371.36 


Si II 


2 


-0.292 


0.069 ± 0.027 


0.130 ± 0.029 


0.116 ± 0.024 


0.137: 


6548.03 


[N II] 


IF 


-0.318 


9.23 ± 0.39 


8.21 ± 0.34 


3.89 ± 0.18 


7.63 ± 0.54 


6562.82 


H I 


Ha 


-0.320 


282 ± 10 


284 ± 10 


288 ± 10 


285.1 ± 9.8 


6578.05 


C II 


2 


-0.322 






0.071 ± 0.019 




6583.41 


[N II] 


IF 


-0.323 


29.2 ± 1.1 


24.61 ± 0.95 


11.83 ± 0.48 


22.5 ± 1.1 


6678.15 


He I 


46 


-0.336 


3.38 ± 0.17 


3.32 ± 0.16 


3.00 ± 0.15 


2.83 ± 0.34 



6 



Table 2 — Continued 



Ao (A) 


Ion 


Mult. 


/(A) 


(A) HII-2 


(B) HII-1 


(C) UV-1 


(D) UV-2 


6716.47 


[S II] 


2F 


-0.342 


13.37 ± 0.55 


11.52 ± 0.47 


13.90 ± 0.57 


36.0 ± 1.5 


6730.85 


[S II] 


2F 


-0.344 


12.44 ± 0.51 


11.03 ± 0.45 


12.17 ± 0.50 


29.0 ± 1.3 


7002.23 


O I 


21 


-0.379 


0.071: 


0.087 ± 0.022 


0.085 ± 0.020 


0.244: 


7065.28 


He I 


1/10 


-0.387 


4.36 ± 0.20 


5.35 ± 0.23 


3.04 ± 0.14 


2.19 ± 0.23 


7135.78 


[Ar III] 


IF 


-0.396 


13.51 it 0.56 


13.07 ± 0.54 


10.40 ± 0.44 


8.89 ± 0.45 


7155.14 


[Fe II] 


14F 


-0.399 


0.047: 


0.053 ± 0.019 


0.048 ± 0.016 


0.094: 


7160.13 


He I 


1/10 


-0.399 




0.015: 






7281.35 


He I 


45 


-0.414 


0.579 ± 0.068 




0.451 ± 0.040 




7318.39 


[O II] 


2F 


-0.418 


2.73 ± 0.15 


2.74 ± 0.13 


2.78 ± 0.13 


4.86 ± 0.32 


7329.66 


[O II] 


2F 


-0.420 


2.09 ± 0.12 


2.10 ± 0.10 


2.23 ± 0.11 


3.56 ± 0.28 


7377.83 


JNi II] 


2F 


-0.425 


0.074: 


0.052 ± 0.019 


0.083 ± 0.020 


0.312: 


7411.61 


[Ni II] 


2F 


-0.429 


0.019: 


0.011: 


0.016: 




7423.64 


N I 


3 


-0.431 


0.011: 


0.010: 


0.007: 




7442.30 


N I 


3 


-0.433 


0.028: 


0.032: 


0.035: 


0.098: 


7452.50 


[Fe II] 


14F 


-0.434 


0.024: 


0.018: 


0.017: 




7468.31 


N I 


3 


-0.436 


0.055: 


0.061 ± 0.020 


0.071 ± 0.019 


0.159: 


7499.85 


He I 


1/8 


-0.439 


0.029: 


0.025: 


0.360 ± 0.036 




7530.60 


C II 


16.08 


-0.443 


0.049: 


0.051 ± 0.018 


0.023: 




7751.10 


[Ar III] 


2F 


-0.467 


3.48 ± 0.18 


3.15 ± 0.15 


2.64 ± 0.13 


2.17 ± 0.23 


8000.08 


[Cr III] 


IF 


-0.492 


0.037: 


0.047 ± 0.017 


0.037: 




8045.63 


[CI IV] 


IF 


-0.497 


0.109 ± 0.042 


0.141 ± 0.025 


0.032: 




8084.00 


He I 


4/18 


-0.500 


0.028: 








8125.31 


Ca I] 


-0.504 


0.012: 








8210.72 


N I 


2 


-0.512 


0.014: 


0.024: 






8216.34 


N I 


2 


-0.513 


0.036: 


0.026: 






8223.14 


N I 


2 


-0.514 


0.046: 


0.041 ± 0.016 






8271.93 


H I 


P33 


-0.518 






0.045 ± 0.015 




8276.31 


H I 


P32 


-0.518 






0.046 ± 0.016 




8281.12 


H I 


P31 


-0.519 






0.059 ± 0.017 




8286.43 


H I 


P30 


-0.519 






0.091 ± 0.020 




8298.83 


H I 


P28 


-0.521 




0.069 ± 0.018 


0.084 ± 0.019 




8306.11 


H I 


P27 


-0.521 




0.046 ± 0.015 


0.066 ± 0.018 




8314.26 


H I 


P26 


-0.522 


0.089: 


0.113 ± 0.022 


0.095 ± 0.020 




8323.42 


H I 


P25 


-0.523 


0.095: 


0.104 ± 0.021 


0.143 ± 0.024 




8333.78 


H I 


P24 


-0.524 




0.692 ± 0.050 


0.903 ± 0.061 




8345.55 


H I 


P23 


-0.525 


0.111 ± 0.042 


0.105 ± 0.021 


0.156 ± 0.025 




8359.00 


H I 


P22 


-0.526 


0.181 ± 0.047 


0.180 ± 0.026 


0.217 ± 0.028 




8361.67 


He I 


1/6 


-0.526 


0.073: 


0.059 ± 0.018 


0.060 ± 0.017 




8374.48 


H I 


P21 


-0.527 


0.190 ± 0.048 


0.182 ± 0.026 


0.214 ± 0.028 




8392.40 


H I 


P20 


-0.529 


0.260 ± 0.052 


0.232 ± 0.029 


0.264 ± 0.031 


0.49 ± 0.14 


8413.32 


H I 


P19 


-0.531 


0.274 ± 0.052 


0.229 ± 0.029 


0.290 ± 0.032 


31.6 ± 1.3 


8437.96 


H I 


P18 


-0.533 


0.319 ± 0.055 


0.300 ± 0.032 


0.320 ± 0.034 


0.137: 


8446.36 


O I 


4 


-0.534 


0.645 ± 0.069 


0.676 ± 0.049 


0.688 ± 0.051 


0.76 ± 0.16 


8467.25 


H I 


P17 


-0.536 


0.373 ± 0.057 


0.375 ± 0.036 


0.369 ± 0.036 


0.35 ± 0.13 


8486.27 


He I 


6/16 


-0.537 


0.019: 


0.018: 






8502.48 


H I 


P16 


-0.539 


0.446 ± 0.061 


0.454 ± 0.039 


0.463 ± 0.041 


0.43 ± 0.14 


8665.02 


H I 


P13 


-0.553 


0.796 ± 0.075 


1.003 ± 0.064 


0.827 ± 0.058 


0.59 ± 0.15 


8680.28 


N I 


1 


-0.554 


0.049: 


0.043 ± 0.016 


0.044 ± 0.015 




8703.25 


N I 


1 


-0.556 




0.024: 


0.016: 




8711.70 


N I 


1 


-0.556 


0.027: 


0.059 ± 0.018 


0.054 ± 0.016 




8718.83 


N I 


1 


-0.557 




0.020: 






8733.43 


He I 


6/12 


-0.558 


0.019: 


0.024: 


0.013: 




8750.47 


H I 


P12 


-0.560 


1.057 ± 0.086 


0.829 ± 0.056 


1.442 ± 0.087 


1.02 ± 0.17 


8845.38 


He I 


6/11 


-0.567 


0.010: 


0.027: 


0.020: 




8862.79 


H I 


Pll 


-0.569 


1.37 ± 0.10 


1.400 ± 0.082 


1.430 ± 0.087 


1.32 ± 0.19 


9014.91 


H I 


PIO 


-0.581 


1.70 ± 0.11 


1.551 ± 0.090 


2.22 ± 0.13 


1.67 ± 0.21 


9068.90 


[S III] 


IF 


-0.585 


25.0 ± 1.3 


25.5 ± 1.3 


21.6 ± 1.1 


17.94 ± 0.82 


9123.60 


[CI II] 


IF 


-0.589 


0.039: 


0.036: 






9229.01 


H I 


P9 


-0.596 


2.54 ± 0.15 


2.63 ± 0.14 


2.64 ± 0.15 


2.65 ± 0.25 
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Table 2 — Continued 



Ao (A) 


Ion 


Mult. 




(A) HII-2 


(B) HII-1 


(C) UV-1 


(D) UV-2 


9530.60 


[S III] 


IF 


-0.618 


71.0 ± 8.5 


63.8 ± 7.0 


51.3 ± 6.2 


43.7 ± 5.2 


9545.97 


H I 


P8 


-0.619 


2.95 ± 0.18 


3.20 ± 0.18 


4.67 ± 0.25 


2.69 ± 0.25 


10031.20 


Ho I 


7/7 


-0.649 


0.214 ± 0.047 








10049.37 


H I 


P7 


-0.650 


6.18 ± 0.35 


6.65 ± 0.36 


5.62 ± 0.31 


11.73 ± 0.60 



E 
o 

CO 
C33 
CD 

o 

X 




3500 
1,0 r 



0,5 



3550 



3600 



3650 



3700 



3750 



3800 



3850 



H9 






He 


[Ne III] 


H5 Hy 




[Olll] HP 


[Olll] 












He 1 [S II] 




















„««^L — 








[Fe III] He 1 [Ar IV] 

1 « ,1 _ 


iw — — ,L~— 

























0,0 

3800 3900 4000 4100 4200 4300 4400 4500 4600 4700 4800 4900 5000 

1,0 r 



0,5 



0,0 





HP 




[0 III] 


He 1 




Ha+ [N II] 


r [sii] 














[0 1] [S III] 


[0 1] 












1 ICI III] 







IL.I. 


He 1 

jL ,j 











4800 
1,0 



5000 5200 5400 5600 



5800 



6000 6200 



6400 



6600 6800 



0,5 



0,0 



[SII] 



He I 



[Arlll] 
[Oil] 



[Arlll] 



[S III] 



6600 6900 7200 7500 7800 8100 8400 8700 9000 9300 9600 9900 

Wavelength (A) 

Fig. 3. — VLT UVES spectra of region B (HII-1). Some of the brightest emission lines are labeled. Telluric lines 
have not been removed and are specially evident in the 6600-9000 A range. Note the gaps in the intervals 5783-5830 
A and 8540-8650 A. 
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He/Hp 
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Fig. 4. — Reddening constant, c{HP), as a function of 
the absorption equivalent widtli in tlie liydrogen lines. 
Each straight line represent the behaviour of the ratio 
of a given Balmer and Paschen line with respect to 

The observed line intensities must be corrected 
for interstellar reddening. This can be done us- 
ing the reddening constant, c{H(3), obtained from 
the intensities of H I lines. However, the fluxes 
of H I lines are also affected by underlying stel- 
lar absorption. Consequently, we have performed 
an iterative procedure to derive both c{H(5) and 
the equivalent widths of the absorption in the hy- 
drogen hues, Wabsi that we use to correct the 
observed line intensities. As we have both H I 
Balmer and Paschen lines, we perform this proce- 
dure separately for them, following the procedure 
explained in Lopez-Sanchez, Esteban & Garci'a- 
Rojas (2006). 

In Figure 4 we plot the reddening constant, 
c{H[3) , versus the absorption equivalent widths in 
the hydrogen hues, Wabs, for the data of region A. 
It can be noted that absorption in the H I Paschen 
lines is practically negligible. This result is also 
found in the rest of analyzed regions. However, 
H I Balmer lines have considerably differences in 
their Wabs , although the c{H(3) derived from them 
is very similar to those derived from Paschen lines. 
We have assumed the average value between the 
c{Hf3) derived from Balmer and Paschen lines as 
representative of the region. This value can be 
used to estimate the appropriate value of Wabs for 
each Balmer line. The c{Hf5) and Wabs that pro- 
vide the best match between the corrected and the 
theoretical Balmer and Paschen line ratios, as well 



as the finally adopted values of c{Hf3) and Wabs 
for the brightest Balmer lines are shown in Table 3. 

4. Emission line profiles 
and gas kinematics 

Our observations confirm the presence of asym- 
metric wings in the line profiles at the center of 
NGC 5253, as were previously reported in echelle 
spectra by Martin & Kennicutt (1995) and Mar- 
lowe et al. (1995). The line-profiles of regions A 
and B show a broader intense component underly- 
ing a narrower one with velocity shifts of the order 
of lower than 10 km s~^. This broad component is 
redshifted with respect to the narrow one in zone 
A, but appears slightly blueshifted in region B. 
Region C seems to show a rather faint broad com- 
ponent. Finally, the line profiles of region D seem 
to be composed by the combination of two sim- 
ilar narrow components with a velocity shift of 
^ 35 km s^^. Consistent values of the velocity 
shift between components were found by Martin 
& Kennicutt (1995) and Marlowe et al. (1995). 

The velocity separation between components 
reaches a maximum at zone D. The maximum 
width of the broad component occurs in regions 
A and B where the full width at half maximum 
(FWHM) reaches values of the order of 110 km 
s~^. Martin & Kennicutt (1995) report consis- 
tent values of the FWHM of this component. The 
velocity shifts and FWHM ratio between broad 
and narrow component for other lines like [O II], 
[O III], [N II], and [S II] are qualitatively simi- 
lar for a given zone. Similar complex emission- 
line profiles and FWHM of the broad components 
were found by Marlowe et al. (1995) and Mendez 
& Esteban (1997) for different samples of dwarf 
starburst galaxies. 

The high signal-to-noise of our spectra has per- 
mitted to perform a spatially resolved analysis of 
the kinematics of the ionized gas observed in the 
center of NGC 5253. It was studied via the analy- 
sis of the centroids of the total line profiles of Ha, 
[O II] AA3726,3729, [N II] A6583 and [S III] A9069 
along the slit position. We have extracted zones 
of 3 pixels wide that corresponds to 0.54" in the 
sky for the Ha and [N II] emission lines and 0.51" 
for the [S III] line. The width of the extracted 
zones for the position of [O II] is 2 pixels, that 
corresponds to 0.50". In Figure 5 we show the 
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Table 3 

h/3 fluxes, reddening coefficients and equivalent widths. 



(A) HII-2 (B) HII-1 (C) UV-1 (D) UV-2 



F(H/3) (xlQi* erg s"! cm-2) 13.45 -|- 

WiUa) (A) 919 

VK(H/3) (A) 234 

wIh-t) (A) 96 

Wins) (A) 49 

iy{He) (A) 27 

VK(H9) (A) 12 

(Ho) (A) 0.0 ± 

(H5) (A) 0.7 ± 

W^af,» (H7) (A) 2.2 ± 

Wabs (He) (A) 1.2 ± 

Wabs (H9) (A) 0.4 ± 

Wabs (Balmcr) (A) 1.3 ± 

c(H/3) (Balmer) 0.22 ± 

Wabs (Paschen) (A) 0.0 ± 

c(H/3) (Paschen) 0.24 ± 

c(H/3) (adopted) 0.23 ± 



0.43 


1 Q c^o _L n /1 
io.OZ ih U.4z 


10. lo ± 0.00 


_l_ r\ nn 
Z.06 ± 0.U9 




1009 


470 


169 




254 


94 


39 




93 


43 


44 




45 


18 


23 




25 


11 


5 




10 


4 


2 


0.2 


0.0 ± 0.2 


0.0 ± 0.2 


0.7 ± 0.2 


0.2 


1.2 ± 0.2 


0.4 ± 0.2 


0.8 dz 0.2 


0.2 


3.9 ± 0.2 


1.1 ± 0.2 


0.8 ± 0.2 


0.2 


1.5 ± 0.2 


0.7 ± 0.2 


0.4 ± 0.3 


0.3 


0.9 ± 0.3 


0.2 ± 0.3 


0.2 ± 0.3 


0.3 


1.7 ± 0.3 


0.8 ± 0.2 


0.6 ± 0.3 


0.02 


0.36 ± 0.03 


0.23 ± 0.03 


0.09 ± 0.02 


0.1 


0.0 ± 0.1 


0.0 ± 0.1 


0.0 ± 0.1 


0.02 


0.39 ± 0.02 


0.27 ± 0.02 


0.10 ± 0.02 


0.02 


0.38 ± 0.03 


0.25 ± 0.03 


0.10 ± 0.02 



position-velocity (PV) diagrams obtained. All the 
velocities are referred to the velocity of the center 
of region B (HII-1). 

Figure 5 indicates that the ionized gas at the 
center of NGC 5253 follows a sinusoidal pattern, 
indicating that its kinematics is not due to pure 
rotation. The behaviour of our PV diagram is in 
agreement (considering the shape and the veloc- 
ity amplitudes) with that obtained by Martin & 
Kennicutt (1995) for their slit positions 2 and 3 
in the central zones of the galaxy. The sinusoidal 
or wave-like pattern may be produced by distor- 
tions produced by the presence of dynamically de- 
coupled gas systems (e.g. Schweizer 1982). This 
explanation would imply that a merging process 
is ongoing in the central zone of the galaxy. An- 
other possibility is that the velocity pattern is pro- 
duced by the outflow from the central starbursts. 
In fact, the wave- like forms produced by super- 
shells are rather common in starbursting dwarf 
galaxies (Martin & Kennicutt 1995; Marlowe et 
al. 1995; Pustilnik et al. 2004). We also find that 
the PV diagrams of the narrow and broad compo- 
nents are similar, it would imply that the outflows 
affect the dynamics of both components in a sim- 
ilar way. It is remarkable that the PV diagrams 
derived from emission lines of different ions are 



10 




-4 I I 

-25 -20 -15 -10 -5 5 10 15 

relative velocity (km s 

Fig. 5. — Position- velocity diagram for the slit po- 
sition observed using Ha, [O II] AA3726,3729, [N II] 
A6583 and [S III] A9069 emission lines. The abscissa 
shows velocities with respect to the center of zone B. 

in excellent agreement except for the [N II] A6583 
line at regions A (HII-2) and B (HII-1). This fact 
could be related to a possible localized nitrogen 
pollution in those zones (we discuss this in §12). 



10 



Table 4 

Physical Conditions of the Ionized Gas 



Diagnostic 



(A) HII-2 



(B) HII-1 



(c 



Te (K) High 



Te (K) Low 



Te (K) 



[Oil] (A3726)/(A3729) 
[SII] (A6716)/(A6731) 
[NI] (A5198)/(A5200) 

[CI III] (A5518)/(A5538) 
[Fe III] 

[Ar IV] (A4711)/(A4740) 

Adopted value'' 

[O III] (A4959+A5007)/(A43f>3) 
[SHI] (A9069+A9532)/(A6312)'' 
[Arlll] (A7136+A7751)/(A5192) 

Adopted value 

[O II] (A3726+A3729)/(A7320+A7330) 
[S II] (A6716+A6731)/(A4069+A4076) 
[Nil] (A6548+A6583)/(A5755) 

Adopted value" 

He I 



660! 
460l 
670 



1 + 130 
-140 
-, + 160 
-210 
+ 550 
<100 

530: 
750 ± 250 

580 ± 110 



11960t^™ 
12330i^3*° 
120001-™ 



600^ 



n+120 

-130 
,.,„ + 170 

g7P + 510 

""J<ioo 
640: 
650 ± 350 
SlOOlf™ 

610 ± 100 



(C) UV-1 
420±100 

,,f)+130 

1140:: 
380: 
350 ± 300 



370 ± 80 



12010; 
11970 
12100 



1+270 
'-300 
1+580 
-'-720 
1+1100 
1900 



10940 



n+240 
-270 

109101=6^° 



(D) UV-2 

270 ± 80 
190+"° 



300 ± 300 



230 ± 70 



109901*^° 
11300+;?™ 



12100 ± 260 12030 ± 260 10810 ± 230 11160 ± 510 



11300 j:™;; 
io88o;?';°o 

110401™° 



n+560 

-680 

n+650 

-720 

11170+**®° 

±±± /U_920 



11330^ 



8680: 



+ 530 
710 



8180 
10410+™° 



600 



10570+^®° 
8380llf3„ 
lOlOOlifoS 



11170 ± 520 11250 ± 490 10530 ± 470 



10270 ± 300 



10200 ± 300 



9000 ± 300 



10350 ± 650 

10600 ± 450 



*^[SIII] A9532 affected by a sky emission line, so we used the theoretical ratio A9532/A9069 = 2.48. 

''ne(ArlV) was not considered in the average, see text. 
'^Te(SII) was not considered in the average, see text. 



5. Physical conditions of the ionized gas 

We have derived the electron temperature. Tg, 
and density, ng, of the ionized gas using several 
emission line ratios. The values obtained for each 
region are compiled in Table 4. All determinations 
were computed with the IRAF task TEMDEN on 
the NEBULAR package (Shaw & Dufour 1995), 
except for the density derived from [Fe III] emis- 
sion lines (see below). We have changed the de- 
fault atomic data of 0+, S+, and S++ included 
in the last version of NEBULAR (February 2004) 
by other datasets that we consider give better re- 
sults. These changes are indicated in Table 4 of 
Garcia- Roj as et al. (2004). 

In order to derive He, we have used several emis- 
sion line ratios between CELs (see Table 4) and the 
intensities of several [Fe III] lines. The [Fe III] den- 
sity was computed from the intensity of the bright- 
est lines of each region (those with errors less than 
30% and not affected by line blending) and apply- 
ing the computations of Rodriguez (2002). All the 
values of rig are consistent within the errors (see 



Table 4) except those obtained from the [Ar IV] 
lines, which are always higher and are expected to 
be representative of the inner zones of the nebu- 
lae. In the case of zones A and B that difference is 
almost one order of magnitude indicating a strong 
density stratification in the ionized gas. 

Once rig was obtained we used it to derive Tg 
using [O III], [S III], [Ar III], [O II], [S II], and [N II] 
line ratios, and we iterated until convergence. We 
have determined the characteristic temperature of 
the helium zone, Te(He I), in the presence of tem- 
perature fluctuations (see §7) using the formula- 
tion of Pcimbert, Peimbcrt, & Luridiana (2002). 
All determinations of Te are included in Table 4. It 
is important to remark that previous Tg determi- 
nations in the literature were restricted to a single 
indicator: Te(0 III). The values of re(0 III) deter- 
mined by Kobulnicky et al. (1997) are consistent 
with those derived from our data. 

We have assumed a two-zone approximation to 
describe the temperature structure of the nebu- 
lae. Therefore, we used the mean of Te(0 III), 
Te(S III) and T'e(Ar III) as the representative tem- 
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pcraturc for high ionization potential ions whereas 
the mean of Te(0 II) and Te(N II) was assumed for 
the fow ionization potential ion temperatures. We 
do not include Te(S II) in the average because its 
value is around 1500 K lower that the ones ob- 
tained from [O II] and [N II] ratios in three of 
the zones. This difference is also reported in the 
Galactic H II region S 311 by Garcia- Rojas et al. 
(2005) and might be produced by the presence of 
a temperature stratification in the outer zones of 
the nebulae. 

6. Ionic Abundances 

The large collection of emission lines measured 
in the spectra has permitted the derivation of 
abundances for an unprecedented large number of 
ions in NGC 5253 from both recombination and 
collisionally excited lines. 

6.1. He+ abundance 

We have measured a large number of He I lines 
in our spectra of NGC 5253, the maximum in re- 
gion A where we detect 34 lines. This permits 
to determine the He"'"/H+ ratio with much better 
precision than in previous studies. We have fol- 
lowed the same method explained in Garcfa-Rojas 
et al. (2005), who applied the maximum likelihood 
method developed by Peimbert, Peimbert, & Ruiz 
(2000) to derive the He+/H+ ratio and some phys- 
ical properties of the ionized gas. In Table 5 wc 
include the He+/H+ ratios for each individual He I 
line and the final adopted average value for each 
region <He+/H+>. In that tabic, wc also include 
the corresponding value of rsssg obtained from the 
maximum likelihood method as well as the corre- 
sponding parameter, which indicates a reason- 
able goodness of the fit obtained by the method in 
the four regions. 

6.2. Ionic abundances from CELs 

The IRAF package NEBULAR has been used 
to derive ionic abundances of 0+, 0++, N+, S+, 
S++, Ne++, Ar++, Ar+3, C1++, and 01+^ for 
each region from the intensity of CELs. The elec- 
tron densities and temperatures for the high and 
low ionization potential ions used are those corre- 
sponding to the aforementioned two-zone scheme. 
The finally adopted ionic abundances are listed 
in Table 6, they correspond to the mean of the 



abundances derived from all the observed individ- 
ual lines of each ion and weighted by their relative 
intensities. Errors were estimated from the uncer- 
tainties in electron density and temperature and 
those associated to the line intensity ratio with 
respect to H/3. Fe++ abundances have been de- 
rived from 6 or 7 [Fe III] lines except in the case 
of zone D, where only three lines were used. The 
lines selected for the derivation of the Fe"'""'"/H+ 
ratios are those not affected by line-blending and 
with line flux uncertainties less than 30%. We 
have used a 34 level model atom that includes 
the collision strengths calculated by Zhang (1996) 
and the transition probabilities given by Quinet 
(1996). Fe++ abundances are also included in Ta- 
ble 6. 

We have detected several [Fe II] lines in our 
spectra, but they are severely affected by fluores- 
cence effects (Rodriguez 1999; Verner et al. 2000). 
Unfortunately, the [Fe II] A8617 line, which is al- 
most insensitive to fluorescence effects, is not ob- 
served because it lies in one of our narrow obser- 
vational gaps. However, we have measured [Fe II] 
A7155, a line which does not seem to be affected 
by fluorescence effects (Rodriguez 1996). We have 
derived Fe+ abundances from this line, assum- 
ing that /(7155)//(8617) ~ 1 (Rodriguez 1996), 
and using the calculations of Bautista & Pradhan 
(1996). The results obtained imply low concen- 
trations of Fe"*" (see Table 6). Due to the faint- 
ness of the [Fe II] A7155 line, and to the assump- 
tion adopted, the derived Fe+ abundances are only 
rough estimates and will not be used in the Fe 
abundance determination. 

We have detected [CI II] lines in two of the ob- 
served zones. However, the C1+/H+ ratio cannot 
be derived from the NEBULAR routines because 
the atomic data of this ion is not included, in- 
stead we have used an old version of the five-level 
atom program of Shaw & Dufour (1995) that is de- 
scribed by De Robertis, Dufour, & Hunt (1987). 
This version uses the atomic data for Cl^ compiled 
by Mendoza (1983), which are rather uncertain 
(Shaw 2003, personal communication). Therefore 
the C1+/H+ ratio given in Table 6 should be inter- 
preted as a rough approximation to the true one. 

6.3. Ionic abundances from RLs 

One of the main goals of this work has been 
the measurement of RLs of heavy element ions in 
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Table 5 
he+/h+ ratios'* 



Line 


(A) Hn-2 


(B) HH-l 


(C) UV-1 


(D) UV-2 


3819.61 


830 ± 166 


704 ± 152 


707 ± 106 




3964.73 


642 ± 88 


603 ± 46 


579 ± 57 




4026.21 


821 ± 43 


775 ± 38 


772 ± 42 


653 ± 108 


4387.93 


676 ± 79 


769 ± 69 


649 ± 81 




4471.09 


876 ± 34 


794 ± 29 


736 ± 27 


676 ± 68 


4713.14 


899 ± 84 


826 ± 61 


691 ± 84 




4921.93 


775 ± 50 


759 ± 43 


722 ± 51 


625 ± 145 


5875.64 


777 ± 29 


840 ± 31 


758 ± 29 


786 ± 48 


6678.15 


836 ± 42 


838 ± 41 


735 ± 37 


741 ± 90 


7065.28 


803 ± 37 


786 ± 34 


729 ± 34 


732 ± 77 


7281.35 


775 ± 91 




632 ± 56 




Adopted'' 


807 ± 16 


791 ± 14 


729 ± 13 


737 ± 34 


•''3889 


9.61 ± 0.74 


11.68 ± 0.70 


7.33 ± 0.70 


1.61 ± 0.92 




18.19 


19.54 


12.67 


3.35 



'''In units of 10 

^The adopted value includes all the relevant uneertainties in emission 
line intensities, ng, and rsggg, and assuming the adopted values. 
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Fig. 6. — Sections of echelle spectra of zones A, B, and C of NGC 5253 showing the hnes of muhiplet 1 of O 11 and 
the WR bump. Note the broad N III A4641 emission hne blended with O II AA4639,4642 emission lines in regions A 
and B. That line is absent in region C. A broad N III A4634 seems to be also present in region A. 
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Table 6 

Ionic abundances from collisionally excited lines 



12 + log(X"/H+) 


(A) HII-2 


(B) HII-1 


(C) UV-1 


(D) UV-2 
4^=0.00 




t^=0.00 t^=0.072±0.027 


4^=0.00 t^=0.050±0.035 


4^=0.00 t^=0.061±0.024 



N+ 


6.62 ± 0.04 


6.81 ± 0.10 


6.55 ± 0.04 


6.67 ± 0.12 


6.30 ± 0.05 


6.48 ± 0.10 


6.61 ± 0.06 


0+ 


7.59 ± 0.06 


7.80 ± 0.11 


7.58 ± 0.05 


7.72 ± 0.12 


7.81 ± 0.06 


8.01 ± 0.11 


8.11 ± 0.08 


0++ 


8.05 ± 0.03 


8.34 


8.07 ± 0.03 


8.26 Ulil 


8.10 ± 0.03 


8.39 ± 0.14 


7.87 ± 0.05 


Ne++ 


7.34 ± 0.05 


7.65 1°:- 


7.36 ± 0.05 


7.56 till 


7.36 ± 0.05 


7 67 +° " 

-0.15 


7.15 ± 0.09 


S+ 


5.70 ± 0.05 


5.89 ± 0.10 


5.64 ± 0.05 


5 7fi +"-12 

0- -0.10 


5.75 ± 0.05 


5.93 ± 0.09 


6.15 ± 0.07 


S++ 


6.45 ± 0.05 


6.77 ± 0.16 


6.44 ± 0.05 


a aA +0.18 
6.64 „o IS 


6.45 ± 0.05 


6 76 +° " 

° -0.14 


6.35 ± 0.07 


C1+ 


4.07: 


4.24: 


4.03: 


4.14: 








C1++ 


4.39 ± 0.08 


4 67 +°-'^'^ 

-0.14 


4.33 ± 0.06 


4 r;i +0.17 

^■^J- -0.14 


4.50 ± 0.07 


4 77 +0-16 


4.19 ± 0.25 


CP+ 


3.73 ± 0.16 


3.96 ± 0.19 


3.85 ± 0.09 


4.00 ± 0.15 


3.30 ± 0.17 


3.53 ± 0.20 




Ar++ 


5.93 ± 0.04 


f- ,o +0.13 

6.18 _(, 11 


5.90 ± 0.04 


6 06 +° " 

O.UD _o.i2 


5.92 ± 0.04 


6.16 ± 0.13 


5.81 ± 0.07 


Ar^+ 


4.88 ± 0.06 


5.18 ± 0.15 


5.00 ± 0.05 


K in +0.17 
-0.15 


4.34 ±0.11 


4 63 +°-^* 




Fe+ 


4.53: 


4.71: 


4.58: 


4.70: 


4.62: 


4.79: 


4.96: 


Fe++ 


5.53 ± 0.08 


5.83 ± 0.16 


5.48 ± 0.08 


5.67 


5.42 ±0.11 


r 71 +0.18 
"■'^ -0.16 


5.89 13 


log(N+/0+) 


-0.97 ± 0.07 


-0.99 ±0.14 


-1.02 ± 0.07 


-1.05 ±0.16 


-1.51 ± 0.07 


-1.53 ±0.14 


-1.51 ± 0.10 



Table 7 

0++ and c++ abundances from recombination lines'* 



(A) HII-2 (B) HII-1 (C) UV 1 

LTE NLTE LTE NLTE LTE NLTE 



O II A46.'i8.86 38: 20: • ■ • ■ • • 47: 24: 

O II A4641.81 23: 27: • ■ • ■ • • 21 ± 7 25 ± 8 

O II A4649.13 12 ± 4 25 ± 8 12 ± 4 25 ± 8 12 ± 4 29 ± 9 

O II A4650.84 34: 16: 36: 17: 42: 18: 

O II A4661.63 21: 12: 20: 12: 47: 26: 

Sum value (all lines) 20 22 17 20 24 25 

0++/H+ adopted value 20 ± 8 18 ± 7 24 ± 10 

12+log(0++/H+) (RLs) 8.30 ± 0.15 8.26 ± 0.15 8.38 ± 0.15 

12+Iog(0++/H+) (CELs) 8.05 ± 0.03 8.07 ± 0.03 8.10 ± 0.03 



C II A4267.15 6 ± 2 5: 6 ± 3 

12+log(C++/H+) (RLs) T-S^to.lt 7.75: T-Slt°o,ll 

12+Iog(C++/H+) (CELs)'' 7Altl \l 7.43 ± 0.17 7.48 ± 0.23 

log(C+ + /0+ + ) (RLs) -0.48j:°-^^ -0.51: -0.57j:g-|^ 

log(C+ + /0+ + ) (CELs)*" -0.62+°;i^ -0.64+° i^ -0.59j:°;^^ 



'^In units of 10 ^. 

•"Value obtained from UV CELs (Kobulnicky et al. 1997). 



NGC 5253, the first dwarf starburst galaxy where 
these kinds of hnes have been unambiguously de- 
tected. Three of the four zones observed show 
the presence of RLs of heavy clement ions in their 
spectra. These lines are those belonging to multi- 
plet 1 of O 11 (see Figure 6) and C II A4267 (see 
Figure 7). All these lines arc produced by pure 
recombination [see Esteban et al. (1998, 2004), 
and references therein] and their intensities de- 
pend weakly on electron temperature and density. 
We have computed the abundances assuming the 
adopted values of Ug and Te(High) given in Ta- 
ble 4 for each zone. The atomic data used and 
the methodology for the derivation of the abun- 
dances from RLs are the same as in Garcia- Roj as 
et al. (2004). The lines of multiplct 1 of O II arc 
not in LTE for densities rie < 10000 cm^"^ (Ruiz 
et al. 2003). We have used the prescriptions given 
by Pcimbcrt, Peimbert, & Ruiz (2005) to calculate 
the appropriate corrections for the abundances ob- 
tained from individual O II lines. These corrected 
abundances show very good agreement with those 
obtained using the sum of the intensities of all 
the lines of the multiplet, which is not affected 
by non LTE effects. Table 7 shows the 0++ /H+ 
and C++/H+ ratios obtained from individual RLs 
as well as the correction for non LTE effects and 
the corresponding sum values for the O II lines. In 



that tabic wc also compare with the same ratios 
obtained from CELs. In the case of the C++/H+ 
ratio, we have compared our determinations based 
on RLs with those obtained by Kobulnicky ct al. 
(1997) from UV CELs for exactly the same zones 
but with a slightly smaller aperture. Consistently 
with the result found for 0++/H+, the C++/H+ 
ratios derived from RLs are systematically larger 
than those derived from CELs. 

7. Abundance Discrepancy 

cind Temperature Fluctuations 

As it can be seen in Table 7, the discrepancy 
between ionic abundances obtained from RLs and 
CELs are between 0.19 and 0.28 dex in the case 
of 0++/H+ and between 0.30 and 0.40 dex in the 
case of C++/H+. Very similar values have been 
obtained for different Galactic and extragalactic 
H II regions. Torres-Peimbert et al. (1980) pro- 
posed that the abundance discrepancy could be 
related to the presence of spatial temperature fluc- 
tuations because of the different functional tem- 
perature dependence of RLs and CELs. Table 8 
shows the different determinations of t"^ and our 
final adopted value. The two upper rows of the ta- 
ble show the parameters that make the 0++ and 
C++ abundances obtained from RLs and CELs to 
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Table 8 

NGC 5253 ^2 PARAMETER 



Method 


(A) HII-2 


(B) HII-1 


(C) UV-1 


0++ (R/C) 
C++ (R/C) 
Adopted 


0.064 ± 0.035 
n nxA +0.027 

0.072 ± 0.027 


0.050 ± 0.035 

0.073: 
0.050 ± 0.035 


0.060 ± 0.030 
062 +0-027 

0.061 ± 0.024 
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Fig. 7. — Section of eclielle spectrum of region A of 
NGC 5253 sliowing the C II A4267 line. 



coincide [denoted as 0++(C/R) and C++(C/R) in 
the table] . We can see that these two values are 
consistent within the errors for all the zones. The 
final adopted value of is the weighted average 
of the individual determinations for each zone. 

8. Total Abundances 

Total abundances have been determined for He, 
C, N, O, Ne, S, CI, Ar, and Fe and are shown in 
Table 9. The non-detection of the nebular He II 
A4686 line in our deep spectra implies the absence 
of a significant amount of He+"'" and O"^"*" in the 
ionized gas of NGC 5253. Then, the use of the 
relation 0/H = 0+/H+ + 0++/H+ is entirely 
justified. For the rest of the elements, we have to 
adopt a set of ionization correction factors (ICFs) 
to correct for the unseen ionization stages. 

For carbon we have adopted the ICF derived 
from photoionization models by Garnett et al. 
(1999). This correction seems to be fairly appro- 
priate considering the high ionization degree of the 
three brightest zones of NGC 5253 where the C II 



line has been detected. 

In the case of neon we have applied the classi- 
cal ICF proposed by Peimbert & Costero (1969), 
that assumes that the ionization structure of 
Ne is similar to that of O. This is a very good 
approximation for high ionization degree ob- 
jects, where a small fraction of Ne+ is expected. 
Martm-Hcrnandez, Schaerer, & Sauvage (2005), 
using A-band NIR spectroscopy, have estimated 
12-hlog(Ne+/H+)=6.46, for a zone located be- 
tween our regions A and B. Assuming this Ne"*" 
abundance for A and B we obtain an ICF of ^1.12, 
somewhat lower than the value of ~ 1.36 that gives 
our adopted ICF scheme. 

We have measured two ionization stages of S, 
and S"*""*", in all regions, but a significant con- 
tribution of S'^'*' is expected. We have adopted 
the ICF given by Stasihska (1978), which is based 
on photoionization models of H II regions and is 
expressed as a function of the 0+/0 ratio. 

We have measured lines of all possible ioniza- 
tion stages of chlorine in two of the zones (A and 
B). In any case, as it can be seen in Table 6, the 
dominant ionization stage is C1++ and the con- 
tribution of CI"*" to the total abundance is rather 
small, at least in the high ionization degree re- 
gions as A, B, and C. Hence, we have neglected 
any contribution of the singly ionized stage for the 
determination of the total abundance of chlorine 
in zone C. In the case of region D, its lower ion- 
ization recommends to take into account the C1+ 
fraction, for that, we have adopted the relation by 
Peimbert & Torres-Peimbert (1977). 

For argon we have determinations of the Ar^'^ 
and Ar'^+ abundances. We find that Ar"*""*" is al- 
most one order of magnitude more abundant than 
Ar^+, indicating that most Ar is in the form of 
Ar"*""*". However, some contribution of Ar+ is ex- 
pected. Martin-Hernandez, Schaerer, & Sauvage 
(2005) have obtained an upper limit for Ar+/Ar++ 



16 



< 0.35 for a zone located between our regions A 
and B. In our case, we have adopted the ICF of 
Izotov, Thuan, & Lipovetski (1994). 

For Fe, we have measured Unes of two stages 
of ionization of iron: Fe"*" and Fe"*""*", but we ex- 
pect an important contribution of Fe^"*". Since 
the Fc+ contribution is small and uncertain, (see 
§ 6.2) the total iron abundance have been obtained 
from the Fe++/H+ ratio and the ICF obtained by 
Rodriguez & Rubin (2005). 

As it was said in the introduction, some previ- 
ous works have found a remarkable nitrogen en- 
richment -and perhaps also helium (Campbell, 
Terlevich, & Melnick 1986)- in some zones of NGC 
5253. Therefore a special care has been taken for 
the determination of the total abundance of these 
two elements. For helium, the aforementioned ab- 
sence of a measurable intensity of the He II A4686 
line indicates that He"*""*" is not important in all 
the zones observed. However, we have to include 
a correction for the presence of neutral helium in- 
side the nebulae. In any case, the high ionization 
degree of regions A, B, and C implies that the 
correction for neutral helium should be small. We 
have used the empirical ICF proposed by Peim- 
bert, Torres-Peimbcrt & Ruiz (1992) for all the 
different zones, although it is expected to be a 
good approximation only for those zones with high 
ionization degree. In the case of zone D, we find 
that the total helium abundance is higher than in 
the other regions. This may be a spurious result 
due to several reasons: a) this zone has compara- 
tively much larger line intensity uncertainties and 
fewer He I lines available, and b) the ionization 
degree of this zone is the lowest of the; obscjrvcid 
regions, and the ICF for helium is comparatively 
more uncertain. Therefore, we will no longer con- 
sider the He/H ratio of region D in our discussion. 
For comparison, we have also used the photoion- 
ization models by Stasiiiska (1990) to estimate the 
ICF for helium in regions A, B, and C, finding that 
the appropriate models for the properties of those 
zones (models C2C1 and C2D1) give very small 
values of the ICF (-1.02). 

In Table 9 we include estimations of the N 
abundance based on three different ICF schemes. 
The first row includes the values assuming the 
standard ICF factor by Peimbert & Costero 
(1969): N/0 = N+/0+, which is a reasonably 
good approximation for the excitation degree of 



the observed regions of NGC 5253. The second 
row includes the total abundances obtained from 
the formulae provided by Mathis & Rosa (1991), 
which give abundances between 0.02 and 0.07 dex 
higher than the ICF of Peimbert & Costero. The 
third estimation has been obtained using the re- 
sults of the photoionization models by Moore, B. 
D., Hester, J. J., & Dufour, R. J. (2004). This last 
determination gives the highest values of N/H, 
about 0.08 and 0.16 dex higher than the ICF of 
Peimbert & Costero. 

In Table 9 we show the total abundances for the 
different zones of NGC 5253 observed for ^^=0.00 
and the t'^ value adopted for each object given in 
Table 8. Our abundances derived from CELs and 
t^=0 are in good agreement with previous determi- 
nations by Walsh & Roy (1989) and Kobulnicky et 
al. (1997). We confirm the remarkable difference 
in the nitrogen abundance between different areas 
of the nebula: regions A and B show higher N/H 
ratios than regions C and D. This difference seems 
to be real because it cannot be accounted for the 
observational uncertainties, it does not depend on 
the ICF scheme used for N, and it is independent 
on the assumption or not of the temperature fluc- 
tuations paradigm. On the other hand, the 0/H 
ratio seems to be slightly higher in zones C and 
D with respect to the other two zones, ~0.10 dex 
in the case of CELs and ~0.13 dex in the case of 
RLs. This difference has not been reported in pre- 
vious works and is of the order of the estimated 
uncertainties. There is also a slight difference of 
the He/H ratio in zones A, B, and C. If this differ- 
ence is real, it would confirm the localized helium 
enrichment found by Campbell et al. (1986) in a 
region encompassing our zones A and B, but with 
much lower values of the abundance difference in 
our case. 

9. Further analysis of the 
emission line profiles 

In §4 we have shown the presence of differ- 
ent velocity components in the emission line pro- 
files of the four zones studied in NGC 5253. We 
have carried out a further analysis of those com- 
ponents, estimating their physical conditions and 
chemical abundances. We have performed a dou- 
ble Gaussian fit to the profiles of the main emis- 
sion lines making use of the Starlink DIPSO soft- 
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Table 9 

Total abundances and abundance ratios 



12 + log(X/H) 




(A) HII-2 




(B) HII-1 




(C) UV-1 


(D) UV-2 


t^=0.00 


t^=n.072±0.027 


(2=0.00 


t2=0.050±0.035 


t^ = o.nn 


t2=0,061±0,024 




Hc^ 


10.97 ± 0.02 


10.95 ± 0.02 


10.95 ± 0.01 


10.94 ± 0.02 


10.93 ± 0.02 


10.91 ± 0.02 


11.08 ± 0.06 


Heb 


10.92 


10.92 


10.91 


10.91 


10.87 


10.87 


10.88 


c 


7.92 ± 0.15 


7.92 ± 0.15 


7.84; 


7.84; 


7.92 ± 0.16 


7.92 ± 0.16 




nC 


7.21 ± 0.07 


7.47 ± 0.17 


7.16 ± 0.07 


7.33 ± 0.18 


6.77 ± 0.08 


7.01 ± 0.17 


6.81 ± 0.11 


Nd 


7.27 ± 0.06 


7.54 ± 0.14 


7.23 ± 0.06 


7.40 ± 0.15 


6.83 ± 0.07 


7.08 ± 0.15 


6.82 ± 0.07 




7.32 


7.61 


7.25 


7.50 


6.93 


7.25 


6.88 


of 


8.18 ± 0.04 


8.45 ± 0.12 


8.19 ± 0.04 


8.37 ± 0.13 


8.28 ± 0.04 


8.54 ± 0.11 


8.31 ± 0.07 


08 


8.42 ± 0.13 


8.42 ± 0.13 


8.37 ± 0.10 


8.37 ± 0.10 


8.53 ± 0.09 


8.53 ± 0.09 




Ne 


7.47 ± 0.06 


7.76 ± 0.20 


7.49 ± 0.06 


7.67 ± 0.22 


7.54 ± 0.07 


7.82 ± 0.20 


7.61 ± 0.09 


S 


6.60 ± 0.07 


6.91 ± 0.15 


6.59 ± 0.05 


6.79 ± 0.13 


6.58 ± 0.05 


6.88 ± 0.12 


6.57 ± 0.07 


Cl'' 


4.62 ± 0.08 


4.86 ± 0.13 


4.59 ± 0.07 


4.75 ± 0.13 








CI' 


4.53 ± 0.07 


4.79 ± 0.10 


4.49 ± 0.06 


4.66 ± 0.12 


4.59 ± 0.08 


4.85 ± 0.14 




Ar 


5.99 ± 0.04 


6.24 ± 0.10 


5.97 ± 0.04 


6.14 ± 0.12 


5.98 ± 0.04 


6.21 ± 0.12 


6.01 ± 0.07 


Fe 


6.08 ± 0.11 


6.43 ± 0.23 


6.01 ± 0.11 


6.23 ± 0.27 


5.82 ± 0.14 


6.17 ± 0.25 


6.06 ± 0.15 


log(X/0) 


c/o 




-0.53 ± 0.18 




-0.53: 




-0.62 ± 0.19 




N/oJ 


-0.91 ± 0.07 


-0.91 ± 0.18 


-1.02 ± 0.07 


-0.97 ± 0.19 


-1.50 ± 0.08 


-1.46 ± 0.18 


-1.49 ± 0.10 


S/O 


-1.58 ± 0.08 


-1.54 ± 0.18 


-1.60 ± 0.08 


-1.58 ± 0.17 


-1.69 ± 0.09 


-1.66 ± 0.16 


-1.74 ± 0.13 


Ne/O 


-0.71 ± 0.08 


-0.69 ± 0.22 


-0.70 ± 0.08 


-0.70 ± 0.24 


-0.74 ± 0.08 


-0.72 ± 0.22 


-0.70 ± 0.15 


Cl/O^ 


-3.65 ± 0.08 


-3.66 ± 0.15 


-3.70 ± 0.08 


-3.71 ± 0.17 


-3.68 ± 0.09 


-3.69 ± 0.17 




Ar/O 


-2.19 ± 0.07 


-2.21 ± 0.15 


-2.21 ± 0.07 


-2.23 ± 0.17 


-2.30 ± 0.08 


-2.33 ± 0.16 


-2.30 ± 0.13 


Fe/O 


-2.10 ± 0.12 


-2.02 ± 0.25 


-2.18 ± 0.11 


-2.14 ± 0.29 


-2.46 ± 0.14 


-2.37 ± 0.26 


-2.25 ± 0.16 



^ICF from Peimbert, Torres-Peimbert &6 Ruiz (1992). 

I— I '^ICF from StasiAska (1990). 

00 

*=ICF from Peimbert &c Costero (1969). 
*^ICF from Mathis & Rosa (1991). 

®ICF from Moore, B. D., Hester, J. J., & Dufour, R. J. (2004). 

^Based on ionic abundances derived from CELs. 

SO++/H+ from RLs and + /H+ from CELs and t^. 

*^C1/H = C1+/H+ + C1 + + /H+ + C1 + 3/H+. 

^ICF from Peimbert & Torres-Peimbert (1977). 

^Assuming the ICF(N) from Mathis & Rosa (1991). 

^Assuming the ICF(Cl) from Peimbert &; Torres-Peimbert (1977). 



Table 10 

Line intensities'' , physical conditions, and abundances'^ of the velocity components 



(A) HII-2 (B) HII-1 (C) UV-1 (D) UV-2 





Narrow 


Broad 


Narrow 


Broad 


Narrow 


Broad 


Narrow 


Broad 


r/~\ TTl OTOl? 

[U iij o7zD 


41 


69 


74 


54 


78 


101 


112 


144 


rO TTl ^79Q 


47 


75 


97 


54 


99 


130 


166 


204 


ro Till 


4 




5 




4 


4 


3 


3 


H/3 4861 


100 


100 


100 


100 


100 


100 


100 


100 


[O III] 5007 


507 


638 


506 


623 


488 


458 


333 


305 


Ha 6563 


282 


282 


282 


282 


282 


282 


282 


282 


[N 11] 6583 


24 


33 


16 


30 


11 


16 


18 


25 


[S 11] 6716 


12 


15 


13 


11 


11 


25 


38 


33 


[S 11] 6731 


11 


14 


11 


11 


10 


21 


27 


29 


C(H/3) 


0.29 


0.18 


0.42 


0.30 


0.37 


0.12 


0.46 


0.05 


rie (cm— 3) 


510 


610 


380 


760 


380 


350 


170 


320 


Te(0 III) (K) 


11800 


12100 


11600 


12300 


10800 


11100 


11200 


10800 


Te(0 11)'= (K) 


11300 


11400 


11100 


11600 


10500 


10800 


10800 


10600 


12 +log(0+ + /H+) 


8.14 


8.09 


8.04 


8.06 


8.13 


8.05 


7.91 


7.92 


log(0+ + /0+) 


0.64 


0.52 


0.36 


0.62 


0.34 


0.19 


-0.01 


-0.15 


log (N+/0 + ) 


-0.87 


-0.90 


-1.31 


-0.83 


-1.53 


-1.44 


-1.47 


-1.45 


12+log(0/H) 


8.23 


8.21 


8.20 


8.15 


8.30 


8.27 


8.22 


8.31 



^Dcrcddcncd line intensity ratios with respect to I(H/5) — 100. 
'^Determined from the intensity of colhsionally excited lines. 
"^Determined from Te(0 III) using the Garnett (1992) calibration. 



ware (Howarth & Murray 1990). The selected 
lines were: [O II] AA3726, 3729, [O III] AA4363, 
5007, H/3, Ha, [N II] AA5755, 6583, and [S II] 
AA6716, 6731. In Table 10 we indicate the emis- 
sion line intensities for the two components in each 
region as well as the reddening coefficient (de- 
rived from the Ha/H/? ratio), the physical con- 
ditions [ue, Te{0 III) and Te{0 II), that was de- 
termined using the Garnett (1992) calibration be- 
tween Te{0 III) and Te{0 II)], some ionic abun- 
dances derived from CELs (0+, 0++ and N+), 
and the 0/H and N+ZO"*" ratios. From the table, 
we can see that most of the values of the different 
parameters are very similar for the narrow and 
broad component of each zone, and also similar to 
the values obtained for the integrated line profile. 
However, we have found some interesting excep- 
tions. First, the reddening coefficient seems to be 
somewhat lower in the case of the broad compo- 
nents in all cases. Secondly, there is a difference 
in the N/0 ratio of the broad and narrow com- 
ponents of zone B. The narrow component has 
a rather normal ratio for a galaxy of its corre- 
sponding O abundance, and the broad one shows 
a higher nitrogen abundance. This would indicate 
that the broad component contains the localized 
N enrichment. However, this behaviour is not ob- 
served in zone A. Here both components show a 
high N/0, similar to that reported for the inte- 
grated profile. 

10. Absorption lines and 
stellar kinematics 

Although our main analysis of the spectra of 
NGC 5253 is focused on emission lines, we report 
the presence of some absorption features. They 
have been identified as Ca II K A3933 and H 
A3968, Mg I A5167, Na I AA5890, 5896 and Ca II 
AA8498, 8662. The strongest absorption lines are 
Ca II H, K, which are unambiguously detected in 
all regions (see Figure 8 for the example of region 
C). We find that the Ca II H, K lines are split- 
ted into several kinematical components (at least 
four), but two of them are due to Galactic ab- 
sorption systems because of their low heliocentric 
radial velocity {vrad ~0 and ~20 km s~^, respec- 
tively). 

The broader absorption component shows a ra- 
dial velocity rather similar to that of the ionized 



NGC 5253- Region C(UV1) 




Ca II K absorptions 



Wavelength (A) 

Fig. 8. — Ca II K A3933 absorptions for region C 
of NGC 5253. Note the existence of several compo- 
nents. The wavelength axis has been corrected for 
the recession velocity of the galaxy. The two Galactic 
components are the bluest ones. 



gas of each region, and can be interpreted as the 
contribution of the stellar populations of NGC 
5253. A second and narrower component with ra- 
dial velocity ~-|-95 km s~^ is also detected in all 
regions. It might be of Galactic origin because 
NGC 5253 is located in the fourth quadrant of 
the Galaxy and radial velocities of the order of 
-1-100 km are common. However, the relatively 
large galactic latitude, b ^ +30°, makes rather 
unlikely its true Galactic nature. Other possibil- 
ities are that this absorption feature is produced 
by a relatively nearby high-velocity cloud or by 
an intergalactic cloud located at ~1.3 Mpc. Using 
the formulae provided by Smoker et al. (2005), we 
find that the column density of this component is 
around 17 times smaller than that associated with 
NGC 5253. 

Data of radial velocities derived from the Ca II 
K absorption feature associated with NGC 5253 
have been included in Figure 9 to compare the 
kinematics of the ionized gas with that of the stel- 
lar component. Notice that Vrad determined from 
Mg I, Na I, and red Ca II absorptions detected 
in C and D are in good agreement with those de- 
rived from Ca II K. Although we only considered 
4 points for constructing the position- velocity dia- 
gram of the absorption features, it is clear that gas 
and stars are kinematically decoupled. The stellar 
component seems to follow a behaviour consistent 
with a rotation pattern. If these velocity differ- 
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Gas from Ha 
Stars from Ca II K 



-30 -25 -20 -16 -10 -5 5 10 16 20 26 30 35 40 45 60 

V (km s"^) 

Fig. 9. — Position- velocity diagram of tlie Ha emis- 
sion line compared to the stellar component kinemat- 
ics derived from the Ca II K absorption features ob- 
served in each region. 



ences between gas and stars are due to the pres- 
ence of a supershell, its velocity relative to the stel- 
lar velocity is '--^40-50 km s^^. This value agrees 
with the expansion velocity of the supershell de- 
tected by Marlowe et al. (1995) in NGC 5253. 

11. The age of the bursts and 

the massive stellar population 

Although previous authors have derived the age 
of the last star-formation burst of NGC 5253, 
we have taken advantage of the high quality of 
our VLT spectra to make a reassessment of it 
through the Ha equivalent width, VF(Ha). We 
have used STARBURST99 (Leitherer et al. 1999) 
spectral synthesis models to estimate the age of 
each region, comparing the predicted Ha equi- 
valent widths with the observed values. We con- 
sider a spectral synthesis model with a metallic- 
ity of Z=0.008, the appropriate for NGC 5253, 
assuming an instantaneous burst with a Salpeter 
IMF, a total stellar mass of 10*^ Mq and a 100 
Mq upper stellar mass. The ages obtained are 
compiled in Table 11, and are in agreement with 
previous results. 

In Figure 6 we can observe the broad emission 
feature of He II A4686 originating from the stellar 
winds of WR stars, which is present in the four 
observed regions of NGC 5253. The He II A4686 
emission line is the most prominent feature of the 
so-called blue WR bump. Unfortunately, the exis- 



0.10 

0.09 

0.08 

0.07 

O 0.06 
+ 

g 0.05 

^ 0.03 
0.02 
0.01 
0.00 



- Z = 0.004 
-Z = 0.008 



W (HP) 

Fig. 10.— WR/(WR-t-0) versus W(H/3) for Schaerer 
& Vacca (1998) models compared with our derived ra- 
tios. 



tence of the red WR bump, that basically consists 
on the C IV A5808 emission line, could not be stud- 
ied because this spectral region is in our observa- 
tional gap between 5783 and 5830 A. In Figure 6 
it is also evident the presence of broad emission 
features of N III A4634 and N III A4640, which 
are characteristic of WNL stars (Smith, Shara, & 
Moffat 1996). 

We have used evolutionary synthesis models by 
Schaerer & Vacca (1998) for O and WR popula- 
tions in young starburst to estimate the number of 
O and WR stars in the regions using the method 
explained in Lopez-Sanchez, Esteban & Rodriguez 
(2004). We have assumed that all the flux of the 
blue WR bump comes from the broad He II A4686 
emission line (we show its extinction-corrected flux 
in Table 11) and a distance of 3.3 Mpc for NGC 
5253. We only need 1 WNL star in regions A 
and B to produce the observed WR spectral fea- 
tures. The number of O stars was derived using 
the procedure given by Guseva, Izotov, & Thuan 
(2000) (their eq. 6). The derived number of O to- 
tal stars and the WR/(WR+0) ratio are quoted 
in Table 11. 

We have found a very good correspondence 
between our results and the predictions from 
Schaerer & Vacca (1998) models. In Figure 10 we 
plot the WR/(WR-hO) ratio versus M^(H/3) for 
the models with Z=0.008 and 0.004. The metal- 
licity of the analyzed regions inside NGC 5253 is 
Z -0.006 for A and B and -0.008 for C and D. 
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Table 11 

Analysis of O and WR populations in observed knots of NGC 5253. 





(A) HII-2 


(B) HII-1 


(C) UV-1 


(D) UV-2 


F (WR Bump) (10-16 erg s"! cm-^) 


11 ± 3 


10 ± 2 


65 ± 18 


1.2 ± 0.7 


L (WR Bump) (103'5 erg s"!) 


2.25 ± 0.66 


2.01 ± 0.46 


13.1 ± 3.6 


0.24 ± 0.14 


WNL eq. stars'* 


1.3 ± 0.4 


1.2 ± 0.3 


8 ± 2 


0.14 ± 0.08 


L (H/3) (10^6 erg s^^) 


271 ± 9 


272 ± 8 


204 ± 7 


51.6 ± 1.9 


07V stars'- 


57 ± 2 


57 ± 2 


43 ± 1 


10.8 ± 0.4 


Adopted Age (Myr) 


3.0 ± 0.1 


2.7 ± 0.1 


4.6 ± 0.1 


5.1 ± 0.1 




~1.0 


~1.1 


~0.4 


~0.3 


O total stars 


59 ± 2 


51 ± 2 


99 ± 5 


36 ± 2 


WR/(WR+0) 


0.023 ± 0.007 


0.023 ± 0.006 


0.07 ± 0.03 


0.004 ± 0.002 



''Assuming that 1 WNL star has a L(He II A4886) = 1.7 X lO'^^' erg s"' (Vacca & Conti 1992). 

''Assuming that 1 07V star has a L(H/3) = 4.76 x 10^"^ erg s"' (Vacca & Conti 1992). 

■^This pajrameter is defined by 77 =07V/0 (Schaerer & Vacca 1998) and depends of the age of the burst. 



12. The localized nitrogen enrichment 

As it was commented in §1, several authors 
have reported localized nitrogen enrichment in ar- 
eas coincident with our regions A and B. Walsh 
& Roy (1989) and Kobulnicky et al. (1997) pro- 
posed that this localized enhancement is due to 
pollution by the winds of massive stars. However, 
for Kobulnicky et al. (1997) a difficulty of this hy- 
pothesis was the apparent lack of accompanying 
He enrichment in the zones where the N/H ratio 
is high. This fact was in contradiction with the ob- 
servational and theoretical understanding of mas- 
sive star yields. In this sense, the slightly higher 
He/H ratio we find in zones A and B with respect 
to C can solve somehow this puzzle. However, it 
must be taken into account that the detection of a 
modest helium pollution is difficult due to several 
reasons: a) the initial abundance of this element 
is far much higher than nitrogen, so the relative 
enrichment should be strong to be unambiguously 
detected; b) the uncertainty introduced by the as- 
sumed ICF to correct for the presence of neutral 
helium. 

Our spectroscopic data include many more He I 
lines than previous works, our measurements have 
much higher signal-to-noisc ratio, and wc arc using 
a more detailed method for determining the He"*" 
abundance. For example, Kobulnicky et al. (1997) 
use a single line: He I A5876, to determine the 
He^/H"*" ratio. The intensity of that line is rather 



affected by collision and radiation transfer effects, 
and the calculations for correcting for these effects 
have improved since Kobulnicky et al. paper. In 
any case, the possible localized helium enrichment 
we detect is only marginal considering the quoted 
uncertainties. Unfortunately, it is clear that any 
result is necessarily unconclusive. 

We have made a rough estimation of the mass of 
newly-created helium and nitrogen (stellar yield) 
necessary to produce the observed overabundances 
in regions A and B. We have reproduced the cal- 
culations of Kobulnicky et al. (1997) for obtaining 
the total ionizing mass of regions A and B, as- 
suming the same angular size and filling factor, 
but considering a distance of 3.3 Mpc instead of 
the 4.1 Mpc assumed by those authors. In Ta- 
ble 12 we show the values we obtain for the cases 
of considering or not temperature fluctuations in 
the abundances, although the difference is almost 
irrelevant. We have considered the mean He and 
N abundances of regions A and B, and the yields 
are computed relative to the abundance of those 
elements in region C. Therefore, we have assumed 
that the initial abundance of NGC 5253 is that 
measured in region C and that zones A and B have 
suffered a localized and very recent increment of 
the Hc/H and N/H ratios. In Table 12 we have 
compared our empirical stellar yields with those 
obtained by Meynet & Maeder (2002) in their stel- 
lar evolution models including rotation. From the 
table, it can be seen that the contribution of a few 
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evolved (WR) massive stars is enough to produce 
the observed pohution, in agreement with the low 
numbers of WNL stars estimated in §11 for zones 
A and B from the flux of the blue WR bump. 
The comparison with the also empirical yields de- 
termined by Esteban et al. (1992) for ring nebu- 
lae associated to Galactic WR stars gives further 
consistency to the hypothesis of the pollution by 
winds of WR stars. In particular, the similarity of 
the ratio of the He and N stellar yields estimated 
for NGC 5253 and for the other objects shown in 
Table 12 is the most remarkable fact of the table. 
That ratio is independent of the quite uncertain 
assumptions considered for deriving the empirical 
yields of the individual elements. 

As it was commented in §8 and can be noted 
in Table 9, the 0/H ratio in zone C is marginally 
higher that in zones A and B, (although the values 
can be considered similar taking into account the 
errors). It is interesting to note that nebular ejecta 
around WR stars (Esteban et al. 1992) and LBV 
stars (Smith et al. 1997; Lamers et al. 2001) show 
a substantial oxygen deficiency in their chemical 
content. This is also consistent with the scenario 
of pollution by evolved massive stars. 

The presence of localized chemical enrichment 
in the youngest starburst of NGC 5253 indicates 
that the timescale of the process should be very 
short, as it was also suggested by Pustilnik et al. 
(2004). The similarity of the enrichment pattern 
observed in this galaxy and that of the Galactic 
WR ring nebulae (with estimated lifetimes of the 
order of 10**-10^ years, see Esteban et al. 1992) 
is another piece of evidence in the same direction. 
The most suitable scenario is that the pollution 
process is produced by the ejection of chemically 
enriched stellar external layers which are pho- 
toionized by the surrounding massive stars. This 
is, for example, the formation mechanism of the 
nebulae around luminous blue variables (LBVs), 
which also show a similar N enhancement (Smith 
et al. 1997; Lamers et al. 2001). In fact, Lamers 
et al. (2001) propose that the LBV nebulae are 
ejected during the blue supergiant phase of the 
progenitor star and that the chemical enrichment 
is due to rotation-induced mixing. Conversely, the 
observed pollution is rather unlikely to be pro- 
duced by stellar wind material confined within a 
hot bubble. In this case, the timescale for mixing 
would be far much larger, of the order of ~100 



Myr (Tenorio-Tagle 1996) because it should cool- 
down to temperatures of the order of the ionized 
gas. 

13. About the origin of the abundance 
discrepancy in NGC 5253 

Very recently, Tsamis & Pequignot (2005) have 
proposed an explanation for the abundance dis- 
crepancy based on the implicit acceptation of 
the existence of temperature fluctuations but pro- 
duced by a chemically inhomogeneous medium. 
These authors propose a chemically inhomoge- 
neous model for 30 Doradus of low temperature 
and high metallicity embedded in a lower den- 
sity medium of high temperature and low metallic- 
ity. The dense high metallicity regions come from 
Type II SN material that has not been mixed with 
the bulk of ISM and is in pressure balance with 
the normal composition ambient gas. These inclu- 
sions would be responsible for most of the emission 
in RLs with virtually no emission in CELs due 
to their very low electron temperature. We will 
present some objections to their models below. 

First, the values of the abundance discrepancy 
factor (ADF) of 0++ (defined as the difference be- 
tween the 0++ abundance determined from RLs 
and CELs) found for very different HII regions in 
different hosts galaxies are rather similar, as it can 
be seen in Table 13. In this table we show the ADF 
values determined for a sample of Galactic and ex- 
tragalactic H II regions where the ADF of 0+"'" has 
been determined, their 0/H ratio, the morphologi- 
cal type and absolute magnitude of the host galaxy 
and the Galactocentric distance in the case of the 
Galactic nebulae. It can be seen that the ADF is 
similar for most of the objects shown in the ta- 
ble, independently of the mass, type and metal- 
licity of their host galaxy. The galaxies included 
in Table 13 have very different metallicities and 
gravitational potential wells and they should also 
have very different star formation histories and SN 
rates. In particular, in the case of NGC 5253, 
a high likelihood of a blow-out has been recently 
postulated by Summers et al. (2004). The energy 
injection rate into this galaxy due to the mechan- 
ical action of massive stars woiild seem siifficicnt 
to allow the expanding hot gas to escape its grav- 
itational potential well. Moreover, the relatively 
small H I halo (~1.8 kpc along the minor axis) of 
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Table 12 

Comparison of empirical and theoretical stellar yields'* 



Sicllar yield 


NGC 5 

(!- = (}) 


253 (A and B) 

('->()) 


Contribution of stellar winds'" 

10 M-:- 60 M :; 
UV,„=0) (r„„=30n) (r„„=0) 


Galactic WR 
King N('l)ulac^ 


mpN 


0.03 


0.04 


0.002 


0.026 


0.015 


0.002-0.02 


"'Pile 


5 


7 


0.2 


1.1 


2.1 


0.2-3.6 


"'/'IK/'"'/'N 




175 


100 


158 


110 


!)n-i8n 



'^Values given in solar masses. 

''Assuming Z / Zq=Q.2 and different initial stellar rotation velocities (in km s^^) (Meynet & 
Maeder 2002). 

"^Empirical yields derived for a sample of Galactic ring nebulae around WN stars (Esteban et al. 
1992). 



Table 13 

Comparison of the abundance discrepancy factor in different HII regions 



Object 


Galaxy 


Galaxy 
Type 


Mv 


Rg" 
(kpc) 


12+log(0/H)b 


ADF(0++)'= 
(dex) 


ADF(C++)'^ 
(dex) 


Reference'' 


M 16 


Milky Way 


Spiral 


-20.9 


6.34 


8.56 


0.45 




1 


M 8 








6.41 


8.51 


0.37 


0.35 


2 


M 17 








6.75 


8.52 


0.27 




2 












8.56 


0.32 




3 


M 20 








7.19 


8.53 


0.33 




1 


NGC 3576 








7.46 


8.56 


0.24 




4 












8.52 


0.32 




3 


Orion nebula 








8.40 


8.51 


0.14 


0.40 


5 












8.52 


0.11 


0.38 


3 


NCC 3603 








8.65 


8.46 


0.29 




1 


S 311 








10.43 


8.39 


0.27 




6 


NCC 5461 


M 101 


Spiral 


-21.6 




8.56 


0.29 


-0.03/0.34 


7 


NCC 604 


M 33 


Spiral 


-18.9 




8.49 


0.20 




7 


30 Doradus 


LMC 


Irregular 


-18.5 




8.33 


0.21 


0.25 


8 












8.34 


0.30-0.43' 


0.41 


3 


NllB 










8.41 


0.69-0.91' 




3 


N66 


SMC 


Irregular 


-16.2 




8.11 


0.36 




3 


Region V 


NGC 6822 


Irregular 


-16.0 




8.08 


0.29 




9 


NGC 2363 


NGC 2366 


Irregular 


-18.2 




7.87 


0.34 


0.31: 


7 


Zones A and B 


NGC 5253 


BCDG 


-17.2 




8.18 


0.23 


0.41S 


10 



*Only for Milky Way objects. 
''Value determined from CELs. 

'^Defined as log(0++/H+)(RLs) - log(0++/H+)(CELs). 
^Defined as log(C++/H+)(RLs) - log(C++/H+)(UV CELs). 

°1- Carci'a-Rojas et al. (2006); 2- Garci'a-Rojas et al. (in preparation); 3- Tsamis et al. (2003); 4- Garci'a-Rojas et al. (2004); 
5- Esteban et al. (2004); 6- Garci'a-Rojas et al. (2005); 7- Esteban et al. (2002); 8- Peimbert (2003); 9- Peimbert et al. (2005); 
10- This work. 

'The highest values of the ADFs include a correction for underlying scattered light. 
^Considering only zone A. 
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this dwarf galaxy would also offer little resistance 
to this blow-out. The likelihood of past blow- 
outs suggests that the proposed delayed enrich- 
ment mechanism based on the rain of coolcd-down 
droplets of SN ejecta would be less efficient in this 
galaxy, but this is not reflected in a lower value of 
the ADFs in NGC 5253 with respect to the other 
H II regions inside more massive galaxies. On the 
other hand, the Galactic H II regions studied in 
our sample arc located at different galactocentric 
distances (from 6 to 11 kpc) in the Galactic disk. 
If the hypothesis of the enrichment by H-poor in- 
clusions of SNe ejecta is correct, this would imply 
that this process is rather independent of the con- 
ditions and properties of the Galactic disk in such 
interval or radial distances. 

Another important piece of evidence in this dis- 
cussion is the fact that the ADF of C++ in NGC 
5253 (and other objects where this quantity has 
been determined, see Table 13) is rather similar 
to that of oxygen. We do not see any metallicity 
dependence of the relative importance of the ADF 
of C++ with respect to that of 0++ in the ob- 
jects where both quantities have been estimated 
(see Table 13). According to the model by Tsamis 
& Pequignot (2005), the overabundances in the 
metal rich regions amount to a factor of 8 in O 
and a factor of 14 in C. It is well known that mas- 
sive star evolution models predict (e.g. Maeder 
1992; Portinari, Chiosi, & Bressan 1998) that most 
of the net C is produced before the onset of the 
SN stage via stellar winds, and since the wind 
strength increases with metallicity, the C yield in- 
creases with metallicity. Carigi, Colin, Peimbert 
(2006) have produced a chemical evolution model 
for NGC 6822, a dwarf irregular galaxy of the Lo- 
cal Group, that predicts that only 36% of the C 
is due to massive stars that produce SN of type 
II, while 63% of the C is due to low and inter- 
mediate mass stars that do not produce Type II 
SN and, therefore, do not contribute to the C con- 
tent of the droplets. Matteucci (2006) also finds 
that for the solar neighborhood most of the C is 
due to low and intermediate mass stars, while Ca- 
rigi et al. (2005) find that about half of the C 
comes from intermediate mass stars and half from 
massive stars. Considering that the 30 Doradus 
chemical composition is similar to that of NGC 
6822 V (see Table 13) one can reasonably infer 
that the expected overabundance of C should be 



lower than that of O, contrary to the assump- 
tions of Tsamis &: Pequignot (2005). Similarly, 
in the chemical evolution models most of the N is 
produced by low and intermediate mass stars and 
not by massive stars, consequently a considerable 
smaller overabundance of N than O is expected, 
contrary to the N and O values adopted in the 
model by Tsamis & Pequignot (2005). 

14. Conclusions 

We present deep echelle spectrophotometry in 
the 3100-10400 A range of four selected zones 
of the nucleus of the blue compact dwarf galaxy 
NGC 5253. We have identified and measured a 
large number of emission lines in the spectra. This 
represents the largest collection of optical emission 
lines available for a dwarf starburst galaxy. Elec- 
tron densities and temperatures have been consis- 
tently determined from a large number of emission 
line ratios of different ions. Chemical abundances 
of He, N, O, Ne, S, CI, Ar, and Fe have been de- 
rived following standard methods. Recombination 
lines of C II and O II have been measured in the 
spectra of three of the individual zones analysed. 
The first time these kinds of lines have been un- 
ambiguously detected in a dwarf H II galaxy. The 
ionic abundances of C++ and 0++ derived from 
these lines are from 0.20 to 0.40 dex larger than 
those determined from the intensity of collision- 
ally excited lines. This behavior has been found 
in all the Galactic and extragalactic H II regions 
where both kinds of lines have been reported for 
the same ion. The quoted abundance differences 
can be accounted for a temperature fluctuations 
parameter {t^) between 0.050 and 0.072. 

The emission line profiles are complex and seem 
to be the combination of —at least— two kine- 
matical components with slightly different radial 
velocity and, in some zones, different line width. 
Apparently, there are no clear differences in the 
physical and chemical properties of the kinemat- 
ical components, except perhaps a somewhat dif- 
ferent reddening coefficient and N/0 ratio in some 
particular zones. Position-velocity diagrams of 
different emission lines show a sinusoidal shape. 
However, the position-velocity diagram of stellar 
absorption spectral features shows a pattern en- 
tirely consistent with rotation. This fact suggests 
that the radial velocity behavor of the ionized gas 



25 



should be due to outflows from the central star- 
bursts and not a product of merging. We have 
derived the ages and the massive star population 
of the starburst. In particular, the detection of 
the so-called blue WR bump indicates the pres- 
ence of a rather small number of WR stars in the 
starbursts. 

Our observations confirm previous results that 
indicate the presence of localized N enrichment in 
two of the central starbursts of NGC 5253. More- 
over, we also detect a possible slight He pollution 
in the same zones. The enrichment pattern is en- 
tirely consistent with that expected by the pollu- 
tion of the winds of masive stars in the WR phase. 
Moreover, the estimated mass of newly created N 
and He is also consistent with the number of WR 
stars determined in the starbursts. 

Finally, we conclude that a recent hypothe- 
sis that tried to solve the abundance discrepancy 
problem in H 11 regions in terms of a delayed 
enrichment by SNe ejecta inclusions (Tsamis & 
Pequignot 2005) seems not to explain the avail- 
able observational data in a satisfactory manner. 

This work is based on observations collected at 
the European Southern Observatory, Chile, pro- 
posal number ESO 70.C-0008(A). We are grate- 
ful to Maria Teresa Ruiz for her collaboration in 
the initial stages of this project. We thank Con- 
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cia y Tecnologfa (MCyT) under project AYA2004- 
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UNAM (grant IN114601). MR acknowledges sup- 
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